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1 Introduction 

Measurements of fermion-pair production in e^e~ collisions at high energies provide a sensitive 
test of Standard Model predictions, and allow limits to be set on many possible new physics 
processes [1-4]. In this paper we present measurements of cross-sections and angular distribu- 
tions for hadronic and Icpton-pair final states at centre-of-mass energies between 189 GeV 
and 209 GeV; forward-backward asymmetries for the leptonic states are also given. The data 
were collected by the OPAL detector at LEP in 1998, 1999 and 2000. 

In the Standard Model, fermion-pair production proceeds via s-channel photon and Z dia- 
grams, except for the e+e" final state where i-channel diagrams dominate. A general feature 
of e"'"e~ collision data at these energies is radiative return to the Z. If one or more initial- 
state radiation photons are emitted which reduce the effective centre-of-mass energy of the 
subsequent e"'"e~ collision, V^, to the region of the Z resonance, the cross-section is greatly 
enhanced. A separation can be made between these radiative events and non-radiative events 
for which x/s' ~ ^/s. While the properties of radiative events are similar to those measured in 
Z decays at LEP 1, modified only by the boost due to recoil against hard initial-state radiation, 
non-radiative events have different properties, refiecting the increased relative importance of 
photon-exchange processes above the Z resonance. At the centre-of-mass energies considered 
here, the contribution of the photon-exchange diagram to the cross-section is about four times 
greater than that of the Z-cxchangc diagram for and r+r^ final states with ~ ^/s. 

The analyses presented here are similar to those already presented at lower energies [1~3]. 
We use identical techniques to measure s' and to separate non-radiative events, which have little 
initial-state radiation, from radiative return to the Z peak. We define non-radiative events as 
those having s' / s > 0.7225, while inclusive measurements correspond to s' / s > 0.01. We 
correct our measurements of hadronic, and r+r^, but not e^e~, events to remove the 

effect of interference between initial- and final-state radiation, as in our previous publications. 
The treatment of the four-fermion contribution to the two-fermion final states is similar to that 
at lower energies. The precise signal definition is discussed in Section 2. 
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While the event selection for hadronic and e"'"e~ final states is essentially unchanged from 
previous analyses, an improvement in the rejection of cosmic ray events in the final 
state has led to a significant reduction in the uncertainty in the residual background. The 
event selection for t~^t~ events has been tightened to reduce the background in this channel, 
also reducing one of the larger systematic uncertainties. In all channels, the higher luminosity 
and hence higher statistics now available have enabled a more thorough study of systematic 
effects. Combining data from three years has led to a significant reduction of the experimental 
systematic uncertainties compared with previous analyses; for example, the systematic error 
on the non-radiative hadronic cross-section has been reduced by ~40%. To take advantage of 
these reduced systematic errors in fits to the Standard Model and searches for new physics we 
present updated results for the high statistics data at 189 GeV, together with new results at 
higher energies; the 189 GeV results supersede those presented in [1]. 

Measurements of fermion-pair production up to 189 GeV have shown very good agreement 
with Standard Model expectations [1-4]. Here we repeat our measurement of the electromag- 
netic coupling constant ae^{^/s) including the higher energy data. In addition, we combine 
results at energies above the Z peak (LEP 2) with those from data taken around the Z peak 
(LEP 1) to determine the mass of the Z boson and the size of the 7-Z interference contribu- 
tion within the framework of the S-matrix formalism [5]. Including data at higher energies 
also allows us to extend the searches for new physics presented in [1]. In particular we obtain 
improved limits on the energy scale of a possible four-fermion contact interaction, and also 
present limits on the mass of a possible heavy Z' boson. 

The paper is organized as follows. In Section 2 we discuss the signal definition, theoretical 
considerations and the corrections made to the data to obtain measurements corresponding 
to this definition. The data and Monte Carlo samples used in the analysis are described in 
Section 3, while Section 4 describes the data analysis and the cross-section and asymmetry 
measurements. In Section 5 we compare our measurements to the predictions of the Standard 
Model and use them to measure the energy dependence of aem- The S-matrix analysis is 
presented in Section 6 and the results of searches for new physics in Section 7. 

2 Signal Definition 

To make precise tests of the Standard Model, the measurements of two-fermion processes must 
be compared with theoretical predictions calculated by, for example, the semi-analytical pro- 
gram Zfitter [6]. We therefore need a signal definition for which the theoretical predictions 
can be made, and also which corresponds closely to the experimental measurements. The defi- 
nition of the two-fermion signal used in this paper is the same as in previous publications [1-3]. 
For the e'^e" final state it is described in Section 2.3 below. For hadronic, /^"'"/x" and r+r" final 
states it is as follows: 

• s' is defined as the square of the mass of the Z/7 propagator. A 'non-radiative' sample 
of events is defined by s'/s > 0.7225, while inclusive measurements correspond to s'/s > 
0.01. 

• Interference between initial- and final-state radiation makes the definition of s' ambiguous. 
To remove this ambiguity, the predicted contribution of interference is subtracted from 
the measured cross-sections, as described in Section 2.1. 
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• Four-fermion final states with a secondary pair arising from an initial-state photon in an 
s-channel diagram (i.e. initial state non-singlet photon diagrams, ISNS^ [7]) are consid- 
ered to be signal if the primary pair passes the s' cut. Four-fermion final states arising 
from final-state photon diagrams (FS-y) are included in the signal. This corresponds to 
definition 1 from [7] page 346. The procedure is discussed in detail in Section 2.2. 

• Cross-section and asymmetry measurements are corrected to full iir acceptance. 

2.1 Interference between Initial- and Final-state Photons 

The data include the effects of interference between initial- and final-state radiation, which 
needs to be subtracted from the measurements to form an unambiguous signal definition for 
comparison with theoretical predictions. We have investigated two methods of performing this 
subtraction. The first method is identical to that used in previous analyses, described fully 
in [3]. We define a differential 'interference cross-section', d^awsR/dm^j dcos6', as the difference 
between the differential cross-section including interference between initial- and final-state radi- 
ation and that excluding interference, as calculated by Zfitter^ . The differential interference 
cross-section may be either positive or negative, depending on the values of the cosine of the 
angle 6 between the fermion and the electron beam direction, and the invariant mass of the 
fermion pair m^j. We then estimate the fraction of this cross-section accepted by our selection 
cuts by assuming that, as a function of cos^ and mfj, its selection efficiency eiFSR(cos^^, mfj) 
is equal to enoint(cos ^, m^j), where enoint has been determined from Monte Carlo events which 
do not include interference. The selected interference cross-section is then subtracted from 
the measured cross-section before efficiency correction, as for other backgrounds. As the ac- 
cepted cross-section is estimated as a function of cos 6', the correction is easily applied to total 
cross-sections, angular distributions or asymmetry measurements. 

The second method investigated uses special samples of Monte Carlo events generated with 
the KK2f program [8] . These samples were generated including initial-final-state photon inter- 
ference with event weights allowing them to be reweighted to exclude the effects of interference. 
Applying selection cuts to these samples allowed the accepted 'interference cross-section' to be 
determined. The accepted cross-sections were in good agreement with those derived from the 
first method. In the case of hadronic events, the average of the two methods was used to correct 
the data and half the difference between them taken as the associated systematic error. For 
li'^li" and t'^t~ final states the first method was used. As at lower energies, the systematic 
error was assessed by repeating the estimate assuming that the efficiency in each bin of cos 6 
and m^j was increased by half its difference with respect to the efficiency at m^f = at the 
same cos^^. 

The corrections to the final measured cross-sections and asymmetry measurements are given 
in Table 1. 

2.2 Four-fermion Effects 

The treatment of the four-fermion contribution to the two-fermion final states is similar to 
that at lower energies. Secondary pairs arising from initial-state photons in s-channel diagrams 

^Cross-sections including interference were calculated by setting the flag INTF=2, those excluding interfer- 
ence by setting INTF=0. 
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(ISNS^) are considered to be signal if the primary pair satisfies the s' cut. Pairs arising from 
final-state photons (FS^) are always considered to be signal. The overall efficiency of event 
selection cuts e is calculated as 



where e^j and t^j^rjr are the efficiencies derived from two-fermion and four-fermion signal Monte 
Carlo events respectively, the generated four-fermion signal cross-section, and cxtot is 

the total cross-section from Zfitter including pair emission. Using this definition of efficiency, 
effects of cuts on soft pair emission in the four-fermion generator are correctly summed with 
vertex corrections involving virtual pairs. For these analyses, a change has been made to the 
method used to separate the signal contribution from the background contribution in the four- 
fermion Monte Carlo events. At lower energies, separate samples of s-channel and t-channel 
four-fermion Monte Carlo events were generated, and the signal contribution was defined by 
kinematic cuts on the s-channel events, designed to include pair production via an initial- 
state photon but to exclude pair production via a Z boson. In this analysis, Monte Carlo 
samples including all four-fermion diagrams, generated with either the grc4f [9] program or 
with KoralW [10] with grc4f matrix elements, were used. Each event was given a weight to 
be signal (or background) calculated using the matrix elements of the appropriate diagrams. 
This method gives a definition of signal which is closer to that employed in the semi-analytic 
calculations with which we compare our results, and also avoids the necessity of generating 
special Monte Carlo samples. 

The inclusion of the four-fermion part of the signal reduces efficiencies by about 0.3% for 
inclusive hadrons, 0.8% for inclusive muons and 1% for inclusive taus. For non-radiative events 
the effects are much smaller, less than 0.02% for hadrons and around 0.2% for muons and taus. 



The discussion above applies to hadronic, /x"*"//" and t+t~ final states. Because of ambiguities 
arising from the t-channel contribution, the acceptance for the e+e~ final state is defined in 
terms of the angle 9 of the electron or positron with respect to the electron beam direction 
and the acoUinearity angle 6'acoi between the electron and positron; a cut on s' is not used. It 
is thus unnecessary to subtract interference between initial- and final-state radiation to make 
an unambiguous signal definition. Cross-sections and asymmetries for e'^e" are not corrected 
for interference between initial- and final-state radiation; they are compared to theoretical 
predictions which include interference. In principle the t-channel process with a second fermion 
pair arising from the conversion of a virtual photon emitted from an initial- or final-state 
electron should be included as signal, as well as the s-channel diagrams. Diagrams with real 
pairs are included in four-fermion Monte Carlo generators, but diagrams with virtual pairs 
are not included, and neither real nor virtual pairs are included in the program we use for 
comparison with the data. It would be improper to treat real pairs alone either as signal or 
background. We choose the best alternative, and simply ignore such events in both efficiency 
and background calculations. If the efficiency for four-fermion events is similar to that for two- 
fermion events, as is expected to be the case here, we are effectively comparing data including 
the four-fermion contribution to theory without. This does not introduce a large error because 
real and virtual pair contributions have opposite sign, and thus their effects tend to cancel in 
any total cross-section. 




(1) 



2.3 e+e- Final States 
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3 Data and Monte Carlo Simulations 



3.1 Data 

The OPAL detector^, trigger and data acquisition system are fully described elsewhere [11-15]. 
The high redundancy of the trigger system leads to negligible trigger inefficiency for all channels 
discussed here. 

The analyses presented in this paper use data recorded during 1998, 1999 and 2000. The 
1998 data were recorded at a centre-of-mass energy near 189 GeV. In 1999, data were taken at 
four different centre- of-mass energy points, close to 192 GcV, 196 GeV, 200 GeV and 202 GeV. 
In 2000, a small amount of data was taken at centre-of-mass energies near 200 GeV and 202 GeV; 
these data have been included with the 200 GeV and 202 GeV data sets taken in 1999. The 
bulk of the 2000 data was taken at a range of centre-of-mass energies between 203 GeV and 
209 GeV, as shown in Fig.l. Also, in 2000, the beam energy was changed within a run by a 
scries of 'miniramps', resulting in a broad distribution of centre-of-mass energies rather than a 
scries of discrete values as in previous years. Data taken while the beam energy was changing 
and other data with a poor beam energy measurement have been removed from these analyses. 
These data amount to about 1.1% of the total 2000 data set. For analysis the good data have 
been divided into two centre-of-mass energy ranges: 202.5 GeV - 205.5 GeV and > 205.5 GeV 
(henceforth referred to as 205 GeV and 207 GeV respectively); these are the energy ranges used 
for the combination of data from all LEP experiments. The mean centre-of-mass energy and 
approximate total integrated luminosity collected at each energy point are shown in Table 2; the 
actual amount of data varies slightly from channel to channel because of differing requirements 
on data quality. 

3.2 Monte Carlo Simulations 

The estimation of efficiencies and background processes makes extensive use of Monte Carlo 
simulations of many different final states. For studies of e+e^ qq we used the KK2f [8] 
program, version 4.13. In KK2f photon radiation is modelled using Coherent Exclusive Ex- 
ponentiation (CEEX) and complete 0{a'^) matrix elements for initial-state radiation are in- 
cluded. Hadronization was performed according to the Pythia6.150 [17] string model. Samples 
hadronized with the Herwig6.2 [18] cluster model or Ariadne4.11 [19] colour dipolc model 
were used for systematic studies. In all cases input parameters have been optimized by a study 
of global event shape variables and particle production rates in Z decay data [20] . Final-state 
radiation from quarks was simulated as part of the hadronization process, and not by the KK2f 
program. For e+e" — > e+e" we used the Bhwide1.04 [21] Monte Carlo program, and for 
e+e"" — > and e+e" t^t~ KK2f was used with KoralZ4.0 [22] for comparison. 

Four-fermion events were modelled with the grc4f [9] generator or with the KoralW [10] 
program with grc4f matrix elements. The latter has superior modelling of initial-state radiation 
in channels without electrons. Final states containing quarks were hadronized using Pythia, 
with Herwig and Ariadne used for systematic studies, as for the e"^e~ — > qq events. Two- 
photon background processes with hadronic final states were simulated using Pythia and 

^OPAL uses a right-handed coordinate system in which the z axis is along the electron beam direction and 
the X axis is horizontal. The polar angle is measured with respect to the z axis and the azimuthal angle (j) 
with respect to the x axis. 
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Phojet [23] at low Q^. At high the TwOGEN [24] program (with the 'perimiss' option [25]), 
Herwig and Phojet were used. In the following, the terms 'tagged' and 'untagged' are used 
to denote the high- and low-Q^ samples respectively. The BDK generator [26] was used to 
simulate two-photon processes resulting in e"'"e~/i+/x~ and e^e~T^r~ final states, while the 
Vermaseren generator [27] was used for the e+e^e+e" final state. The e+e~ 77 background 
in the e+e" final state was modelled with the Radcor [28] program, while the contribution 
from e"'"e~7 where the photon and one of the charged particles are inside the detector acceptance 
was modelled with Teegg [29]. 

Monte Carlo samples were generated at 189 GeV, at the four nominal energy values of the 
data collected in 1999, and at several energies spanning the range 204 GeV to 208 GeV for 
simulation of the data taken in 2000. All samples were processed through the OPAL detector 
simulation program [30] and reconstructed in the same way as for real data. Efficiencies and 
backgrounds at the centre-of-mass energy values corresponding to the data were determined by 
fitting the energy dependence of these quantities. 

For the measurement of the luminosity, the cross-section for small-angle Bhabha scattering 
was calculated using the Monte Carlo program Bhlumi [31], using generated events processed 
through a program which parameterizes the response of the luminometer [15]. 

4 Cross-section and Asymmetry Measurements 
4.1 Luminosity 

The integrated luminosity was measured using small-angle Bhabha scattering events, e+e^ — > 
e'^'e", recorded in the silicon-tungsten luminometer [15]. The luminometer consisted of two 
finely segmented silicon-tungsten calorimeters placed around the beam pipe, symmetrically on 
the left and right sides of the OPAL detector, 2.5 m from the interaction point. Each calorimeter 
covered angles from the beam between 25 mrad and 59 mrad. The luminosity determination 
closely followed the procedure used for the precise determination at LEP 1 [15]. However, before 
LEP 2 data-taking, tungsten shields designed to protect the tracking detectors from synchrotron 
radiation were installed. These introduced about 50 radiation lengths of material in front of 
the calorimeter between 26 mrad and 33 mrad from the beam axis, thus reducing the useful 
acceptance of the detector. 

Bhabha scattering events were selected by requiring a high energy cluster in each side of 
the detector, using asymmetric acceptance cuts. The energy in each calorimeter had to be at 
least half the beam energy, and the total energy in the fiducial region of both calorimeters 
had to be at least three quarters of the centre-of-mass energy. The two highest energy clusters 
were required to be back-to-back in 0, | |0r — 0z,| — 7r| < 200 mrad, where 0r and 0l are the 
azimuthal angles of the cluster in the right- and left-hand calorimeter respectively. They were 
also required to be collinear, by placing a cut on the difference between the radial positions, 
Ai? = Rr—Rl, at \AR\ < 2.5 cm, where Rr and Rl are the radial coordinates of the clusters on 
a plane approximately 7 radiation lengths into the calorimeter at 2; = ±246.0225 cm. This cut, 
corresponding to an acoUinearity angle of about 10.4 mrad, effectively defines the acceptance for 
single-photon radiative events, thus reducing the sensitivity of the measurement to the detailed 
energy response of the calorimeter. The distribution of AR for the data taken in the year 2000 
is shown in Fig. 2(a). 
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Inner and outer radial acceptance cuts delimited a region between 38 mrad and 52 mrad on 
one side of the calorimeter, while for the opposite calorimeter a wider zone between 34 mrad and 
56 mrad was used. Two luminosity measurements were formed with the narrower acceptance 
on one side or the other side. The final measurement was the average of the two and has no 
first order dependence on beam offsets or tilts. The distributions of the radial coordinates of 
the clusters for the data taken in the year 2000 are shown in Fig. 2(b,c). 

The acceptance A of the luminosity measurement is affected by any change in the inner and 
outer edges of the acceptance as follows: 

M « X 10- (2) 

and 

^ . +^ . 10-, (3) 

where Rm and i?out denote the radial coordinates of the inner and outer cuts. The coefficients 
in the expressions given above are determined by simple analytic calculations, using the 1/6^ 
Bhabha spectrum, the nominal half distance between the reference planes of the two calorime- 
ters and the inner and outer acceptance radii (9.45 cm and 12.7 cm). The residual bias on 
the inner and outer cut positions was estimated as at LEP 1, by a procedure called anchoring, 
which is fully explained in [15]. In this approach the fundamental tool is the radial position of 
the silicon pad with maximum signal in a given longitudinal layer. As the radial position of the 
incoming particles crosses a radial pad boundary in a single layer, the average pad-maximum 
moves rapidly from one pad to the next, giving an image of the pad boundary, as shown in 
Fig. 3. The coordinate offset at the inner cut position was determined to be about 30 //m in 
both the right and the left calorimeter; this offset can be seen in Fig. 3. By applying Equa- 
tion (2) this is equivalent to an acceptance variation of +0.14 %. At LEP 1 the absolute value 
of the coordinate offset at the inner cut was less than 10 /im. The larger value in LEP 2 data is 
attributed to the effects of about 2 radiation lengths of prcshowering material, consisting of ca- 
bles and beam pipe support structures, in front of the central angular region of the calorimeter 
including the position of the inner radial cut. The estimated systematic correction was applied 
to the measurement, but the full size of the effect was conservatively kept as a systematic error. 

The errors on the luminosity measurement at each energy are summarized in Table 3. The 
experimental systematic error is dominated by the uncertainty of the inner radial cut. Among 
the other errors are trigger efficiency (0.06%), energy response particularly in the low energy tail 
(0.03%), beam parameters (0.02%), backgrounds (0.02%) and Monte Carlo statistics (0.08%), 
where all the numbers refer to year 2000 data but are fairly similar in the other data samples. 
The error on the theoretical prediction of the Bhabha cross-section of 0.12% is taken from [32]. 

The errors on luminosity are included in the systematic errors on cross-section measurements 
presented in this paper. Correlations between cross-section measurements arising from common 
errors in the luminosity have been taken into account in the interpretation of the results. 
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4.2 Hadronic Events 



4.2.1 Event Selection 

The selection of hadronic events is identical to previous analyses [1-3]. For both inclusive and 
non-radiative samples, the selection efficiency is typically ~85% and the purity is ~92%. 

• To reject leptonic final states, events were required to have high multiplicity: at least 
7 electromagnetic clusters and at least 5 tracks satisfying standard quality criteria [3]. 

• Background from two-photon events was reduced by requiring a total energy deposited 
in the electromagnetic calorimeter of at least 14% of the centre-of-mass energy: i?vis = 
TjEcIus/^/s > 0.14, where Edus is the energy of each cluster. 

• Any remaining background from beam-gas and beam-wall interactions was removed, and 
two-photon events further reduced, by requiring an energy balance along the beam di- 
rection which satisfied i?bai =| ^(ii^cius " cos6'cius) I /^^-^'cius < 0.75, where 6'cius is the polar 
angle of the cluster. 

• At these centre-of-mass energies, the cross-section for production of W"'"W~ events is 
comparable to that for non-radiative qq events, and the above selection cuts have high 
efficiency for those W+W" events with hadrons in the final state. Events selected as 
W"'"W~ candidates using the criteria described in [33] (with reference histograms updated 
for the higher energy data) were therefore rejected. This cut also removes some of the 
(much smaller) contribution from ZZ final states. As a cross-check, an analysis was also 
performed in which W+W" candidates were not rejected, but their expected contribution 
subtracted. 

• The effective centre-of-mass energy, x/s', of the e+e^ collision, determined as described 
below, was required to satisfy s'/s > 0.01 for the inclusive sample and s'/s > 0.7225 for 
the non-radiative sample. 

Distributions of selection variables, before applying the W"'"W" rejection and s' cuts, are shown 
in Fig. 4. The Monte Carlo modelling of these variables is generally very good, except at very 
low multiplicities. 

The effective centre-of-mass energy x/s' of the e+e~ collision was estimated as follows. The 
method is the same as that used in previous analyses [1-3]. Isolated photons in the electro- 
magnetic calorimeter, with a minimum energy of 10 GeV, were identified, and the remaining 
tracks, electromagnetic and hadron calorimeter clusters formed into jets using the Durham {kx) 
scheme [34] with a jet resolution parameter ^cut = 0.02. If more than four jets were found the 
number was forced to be four by adjusting the jet resolution parameter. The jet energies and 
angles were corrected for double counting of energy using the algorithm described in [35] . The 
jets and observed photons were then subjected to a series of kinematic fits imposing the con- 
straints of energy and momentum conservation, in which zero, one, or two additional photons 
emitted close to the beam direction were allowed. The fit with the lowest number of extra 
photons which gave an acceptable was chosen. The value of was then computed from 
the fitted four-momenta of the jets, i.e. excluding photons identified in the detector or those 
close to the beam direction resulting from the fit, which were assumed to arise from initial-state 
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radiation. If none of the kinematic fits gave an acceptable x^, \fs' was estimated directly from 
the angles of the jets as in [36] . The distribution of for all energies combined is shown in 
Fig. 9(a). 

The efficiency of the selection cuts was determined from Monte Carlo events generated with 
the KK2f program, without inclusion of interference between initial- and final-state photon 
radiation, and corrected for the effect of the four-fermion signal component as described in 
Section 2.2. The feedthrough of events with lower s' into the non-radiative sample and expected 
backgrounds were also determined from Monte Carlo. For both efficiencies and backgrounds 
a linear fit to values at ccntrc-of-mass energies between 189 GcV and 208 GeV was used to 
determine the value at the mean centre-of-mass energy of the data. The KoralW Monte Carlo 
events used to estimate the four-fermion background do not include complete electroweak 0{a) 
corrections to the e+e~ — > W+ W~ process which are now available in the program KandY [37] . 
The backgrounds calculated using KoralW were corrected using samples of events generated 
with KandY. These corrections are small for the standard analysis, where W^W~ candidates 
are rejected, but significant for the cross-check analysis in which W+W" events are not rejected. 
In this case they increase the measured inclusive and non-radiative cross-sections by 0.5% and 
1.2% respectively. Efficiencies and backgrounds are summarized in Table 4. 

To measure the angular distribution of the primary quark in the hadronic events, we have 
used as an estimator the thrust axis for each event determined from the observed tracks and 
clusters. The angular distribution of the thrust axis was then corrected to the primary quark 
level using bin-by-bin corrections determined from Monte Carlo events. For the bin size chosen, 
the bin-to-bin migration of events is ~10% to either side. No attempt was made to identify the 
charge in hadronic events, and thus we measured the folded angular distribution. 

4.2.2 Systematic Uncertainties 

The systematic errors on the hadronic cross-sections are summarized in Table 11, with a detailed 
breakdown at 200 GeV given in Table 12. Where no dependence on energy or year was expected 
or seen, the values were determined by combining data at all energies. The resulting high 
statistics have resulted in a reduction of many contributions compared with previous analyses. 

Initial-state radiation modelling. Efficiencies were calculated using the KK2f generator 
with 0{a^) Coherent Exclusive Exponentiation (CEEX) of radiation. To assess the effect 
of initial-state radiation on the selection efficiencies and s' determination, the events were 
reweighted to 0{q) CEEX. In accordance with the recommendations of [8], half the difference 
between 0{a) and 0{a^) was assigned as the systematic error, reflecting the effects of missing 
higher order terms in the perturbative expansion. 

Fragmentation modelling. The effect of the hadronization model on the selection ef- 
ficiencies has been investigated by comparing the string model implemented in Pythia with 
the cluster model of Herwig [18] and the colour-dipole model of Ariadne [19]. To reduce 
the statistical errors on this comparison, the same primary quarks generated with KK2f were 
fragmented according to each model in turn, and the selection efficiencies compared. The de- 
viations of the two predictions from the Pythia value were evaluated. Statistically significant 
differences were seen, and the larger of these was assigned as the systematic error. In addition, 
the effects on the efficiencies of changing the cuts on the number of tracks and clusters by one 
unit were also taken into account, to cover imperfections in the modelling of low multiplicity 
jets. 
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Detector effects. The selection of inclusive events is mainly based on the electromag- 
netic calorimeter, and is thus sensitive to the energy scale of the calorimeter, and any angular 
dependence of the energy scale. For non-radiative events, the selection is sensitive to jet and 
photon energies, angles, and their errors, and jet masses, which are used as input to the kine- 
matic fits used to determine s'. Studies of calibration data taken at the Z peak have been 
used to determine small year-dependent corrections to these parameters in the Monte Carlo 
simulations. Variations of these corrections by their errors were used to assign corresponding 
systematic errors on the cross-sections. The uncertainty in the energy scale of the electromag- 
netic calorimeter leads to an error on the inclusive cross-sections of about 0.2%. The largest 
effect on the non-radiative cross-sections (0.12%) arises from the jet energy scale. 

s' determination. Possible systematic effects in the determination of s' not already covered 
by the studies of initial-state radiation modelling, fragmentation modelling and detector effects 
were studied using two alternative methods of calculating s'. Firstly the default algorithm 
was modified to allow only a single radiated photon, either in the electromagnetic calorimeter 
or along the beam axis. Alternatively the cuts defining photon candidates in the detector 
were varied. The differences, averaged over all centre-of-mass energies, were not statistically 
significant, and the statistical precision of this test was included as a systematic error associated 
with the s' determination. 

W+W~ rejection cuts. The effect of the W+W~ rejection cuts on the signal efficiency 
was studied using distributions of variables which distinguish W^W~ events from qq events. In 
the case of W+W~ qqqq events the QCD matrix element for four-jet production, W420 [38] 
was used. This is an event weight formed from the 0(0.^) matrix elements for the four-jet 
production processes e"'"e~ ^ qq ^ qqqq, qqgg- The distribution of VF420 after all event 
selection cuts except the W+W~ veto is shown in Fig. 5(a) for non-radiative events from 
the combined data sample. A clear separation between signal events and background events 
is seen. Figures 5(b) and (c) show the distributions for those events rejected by and passing 
the W"'"W~ veto respectively. The distribution obtained from the data, before applying the 
W+W" rejection, was fitted with the sum of expected signal and background contributions, 
allowing the absolute normalization of both to vary. The fit region was chosen to include the 
majority of signal events rejected by the W+W~ veto and a majority of background events 
which pass the veto, as illustrated in Fig. 5. The resulting scale factor for the signal, or its 
statistical error, was applied to the rejected signal cross-section to estimate the corresponding 
uncertainty in efficiency. In the case of W+W^ — > qq£z/ events, a similar procedure was applied 
to the distribution of the magnitude of the vector sum of transverse momenta for all visible 
particles. In all cases, the scale factors were found to be consistent with unity. 

Backgrounds. The uncertainty in the W+W" background was estimated from the fits to 

W420 and transverse momentum distributions (for W^W^ qqqq and W"'"W~ qqiu respec- 
tively) described above for the W^W~ rejection cuts. In addition, the effects of initial-state 
radiation modelling and fragmentation modelling on the four-fermion background were inves- 
tigated by varying these models. In the inclusive sample, the largest background uncertainty 
arises from tagged two-photon events. This was investigated by comparing the predictions of 
different Monte Carlo generators: either a combination of Herwig for single-tagged events plus 
Phojet for double-tagged events, or TWOGEN for both. An average of the two prescriptions 
was found to give the best representation of the data at low s', and was used in the cross-section 
determination, with half the difference between the two predictions taken as the systematic er- 
ror. The (small) differences between the Pythia and Phojet programs were used to assess 
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the systematic uncertainty in the untagged two-photon background. Similarly, the small un- 
certainties in the t+t~ background were estimated by comparing the predictions of KK2f and 
KoralZ. 

Interference. The error arising from the subtraction of interference between initial- and 
final-state photon radiation was estimated from the difference between the two methods of 
determination as described in Section 2.1. 

The results of the cross-check analysis, in which events identified as W+W~ were not re- 
moved, but the expected contribution subtracted, are in excellent agreement with those of the 
primary analysis, with slightly larger total errors. 

4.3 Muon Pairs 
4.3.1 Event Selection 

The selection of events is essentially the same as in previous analyses [1], except that a 

small improvement has been made in the rejection of cosmic ray events leading to a reduction in 
the uncertainty associated with this background. The efficiency of the selection cuts is typically 
~74% for inclusive events and ~88% for non-radiative events. The corresponding purities of 
the selected samples are ~90% and ~97% respectively. 

• Muon pair events were required to have at least two tracks with momentum greater than 
6 GeV, I cos ^1 < 0.95, separated in azimuthal angle by more than 320 mrad, and identified 
as muons. These tracks must have at least 20 hits in the central tracking chambers and 
the point of closest approach to the nominal beam axis must lie less than 1 cm in the r-0 
plane and less than 50 cm along the beam axis from the nominal interaction point. To 
be identified as a muon, a track had to satisfy any of the following conditions: 

— At least 2 muon chamber hits associated with the track within an azimuthal angular 
range A(f) = (100 + 100/p) mrad, with the momentum p in GeV; 

— At least 4 hadron calorimeter strips associated with the track within an azimuthal 
angular range A(j) = (20 -|- lOO/p) mrad, with p in GeV. The average number of strips 
per layer, taken over all layers with at least one hit, had to be less than 2 to discrim- 
inate against hadrons. For | cos^| < 0.65, where tracks traverse all 9 layers of strips 
in the barrel calorimeter, a hit in one of the last 3 layers of strips was required; 

— Momentum p > 15 GeV and less than 3 GeV electromagnetic energy associated to 
the track within a cone of half-angle 200 mrad. 

If more than one pair of tracks satisfied the above conditions, the pair with the largest 
scalar sum of momenta was chosen. No requirement was made that the tracks have 
opposite charge. 

• Background from high multiplicity events was rejected by requiring that there be no other 
track in the event with a transverse momentum (relative to the beam axis) greater than 
1.65% of the beam energy. 
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• Background from cosmic ray events was removed using the time-of-flight (TOF) counters 

and vertex cuts. Figure 6(a) shows the distribution of the time difference, At, between 
pairs of back-to-back TOF counters for fi^fi^ candidates, clearly showing one peak at the 
origin from muon pairs and a second peak at about 15 ns from cosmic rays. Events were 
accepted if they had — 20 ns < At < 8 ns and at least one of the time measurements was 
within 10 ns of that expected for a particle coming from the interaction point. If only 
one TOF hit was recorded, it had to be within 10 ns of the expected time. In addition, 
events were required to pass loose cuts on the matching of the central detector tracks to 
the interaction vertex. Events without a good TOF hit were required to pass tight vertex 
criteria. 

• Background from two-photon events was rejected by placing a cut on the total visible 
energy, Ey^s, defined as the scalar sum of the momenta of the two muons plus the energy 
of the highest energy cluster in the electromagnetic calorimeter: 

Ryi, = £;vis/Vs > 0.5(m|/s) + 0.35. 

The value of this cut is 0.15 below the expected value of it!vis for muon pairs in radiative 
return events where the photon escapes detection, visible as a secondary peak in Fig. 6(b). 
Furthermore, for inclusive events, if the ratio of the visible energy to the centre-of-mass 
energy was less than 0.5(m|/s) -|- 0.75 the muon pair invariant mass was required to 
be greater than 70 GeV. For all non-radia tive events the muon pair invariant mass was 
required to be greater than ^/(ml -|- 0.1s). 

• The effective centre-of-mass energy of the e^e~ collision, determined as described 
below, was required to satisfy s' / s > 0.01 for the inclusive sample and s'/s > 0.7225 for 
the non-radiative sample. 

Roughly 10% of selected events have a photon detected in the electromagnetic calorimeter 
with an energy above 30 GeV, separated from the nearest muon by at least 20°. If such an event 
was planar, i.e. the sum of the angles between the three particles (two muons plus photon) was 
greater than 358°, the photon was assumed to be initial-state radiation and s' was calculated 
from the angles of the two muons and the photon using three-body kinematics. For all other 
events the value of s' was estimated from the polar angles 9i and 62 of the two muons, assuming 
massless three-body kinematics to calculate the energy of a possible undetected initial-state 
photon along the beam direction as 

E^^y/^-\ sin(^i + ^2) 1/(1 sm{0i + ^2) | + sin 0^ + sin ^2) ■ (4) 

The observed distribution of ^/s' for all data combined is shown in Fig. 9(b). 

The selection efficiencies and feedthrough of events from lower s' into the non-radiative sam- 
ples were determined from Monte Carlo events generated with KK2f without interference be- 
tween initial- and final-state radiation, corrected for the four-fermion contribution as discussed 
in Section 2.2. Backgrounds were also determined from Monte Carlo simulations. Efficiencies 
and backgrounds at each energy are summarized in Table 4. 

In approximately 2% of events the two muon tracks have the same charge; for the 

asymmetry and angular distribution measurements this charge ambiguity was resolved using the 
acoplanarity of track segments reconstructed in the muon chambers. Acoplanarity is defined as 
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101 — 02 1 — 180° where 0i and 02 are the azimuthal angles of the muon segments. Bending of the 
charged particle trajectories in the magnetic field results in positive or negative acoplanarity 
depending on the charge of the particle with the lower value of 0. In the measurement of both 
the angular distributions and asymmetries, the final values were obtained by averaging the 
distribution measured using the negative muon with that using the positive muon; although this 
averaging does not reduce the statistical errors on the measurements, it is expected to reduce 
most systematic effects. The forward-backward asymmetries at each energy were obtained by 
counting the numbers of events in the forward and backward hemispheres, after correcting 
for background and efficiency. The asymmetries were corrected to the full angular range by 
applying a multiplicative correction obtained from Zfitter to the asymmetry measured within 
the acceptance of the selection cuts (| cos^| < 0.95). 

4.3.2 Systematic Uncertainties 

Systematic errors on the ^~^fi~ cross-sections are summarized in Table 11, with a detailed 
breakdown at 200 GeV given in Table 13. The main contributions arc discussed below. 

Efficiency. The systematic uncertainty in the efficiency was evaluated using high statistics 
LEP 1 data and Monte Carlo samples. The cross-section at the Z peak is well known: it 

has been measured with a systematic uncertainty of about 0.2% [39]. The selection cuts 

were applied to the LEP 1 data and Monte Carlo samples. A statistically significant difference 
between the number of data events selected and the number expected from Monte Carlo was 
observed, and this difference was used to estimate the systematic error associated with the 
efficiency. Most kinematic cuts are a function of ^/s and scale smoothly to the Z peak; for this 
comparison it was necessary only to relax the cut on the visible energy so that the efficiency for 
events on the Z peak remained high. LEP 2 events have a different angular distribution from 
LEP 1 events, and in particular radiative events are boosted towards the endcap regions of the 
detector. The agreement between data and Monte Carlo was therefore checked as a function 
of cos 6*, and the results were reweighted to the angular distribution of the high energy data 
to obtain the systematic error on the total cross-section. To check for possible changes of the 
detector response with time, this procedure was repeated with the calibration data taken at 
the Z during 1998-2000; the observed difference between data and Monte Carlo was consistent 
with that determined from the LEP 1 study, but with poorer statistical precision. 

Initial-state radiation modelling. The systematic error on efficiency derived from the 
LEP 1 data does not include the effects of uncertainties in the modelling of initial-state radia- 
tion. As for hadronic events, this uncertainty was estimated by reweighting KK2f events from 
0{a'^) to 0{a) CEEX and taking half the predicted change in efficiency. 

Feedthrough. The uncertainty in the feedthrough of events with lower s' into the s'/s > 
0.7225 sample was estimated by comparing the prediction of KK2f with that of KoralZ. 

Cosmic background. The uncertainty due to any remaining cosmic background in the 
muon pairs was estimated from the vertex distribution of events after relaxing some of the 
time-of-flight and vertex criteria. 

Other backgrounds. The main backgrounds in the muon pairs arise from various leptonic 

four-fermion final states and from tau pairs. The four-fermion backgrounds are principally 
channels including at least two muons, and include a significant contribution from production 
of e^e^/x+/i^ final states via two-photon processes. Backgrounds were studied by considering 
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distributions of selection variables after loosening some of the selection cuts. The numbers of 
events in data and Monte Carlo were compared for a region enriched in a particular background, 
and the difference, or its statistical error, whichever was greater, used to estimate the systematic 
error from that background source. For example, the two-photon e+e~/x"'"/x~ background was 
studied using the distribution of visible energy after removing the cuts on visible energy and 
muon-pair mass, shown in Fig. 6(b); the comparison was made in the visible energy range 
between 10% and 40% of the centre-of-mass energy, which is completely dominated by this 
background. 

Interference. The uncertainty arising from the removal of the contribution from interfer- 
ence between initial- and final-state radiation was estimated as described in Section 2.1. 

Asymmetry. Systematic uncertainties in the asymmetry measurement were assessed by 
comparing results obtained using different combinations of tracking and muon chambers to 
measure the muon angles. The change in asymmetry when same-sign events were excluded 
from the sample was included as a systematic error. Other small contributions arise from 
the efficiency and background correction and subtraction of interference between initial- and 
final-state radiation. 

4.4 Tau Pairs 
4.4.1 Event Selection 

The selection of e"'"e~ t^t~ events is based on that used in previous analyses [1], using 
information from the central tracking detectors and electromagnetic calorimetry to identify 
events with two coUimated, low multiplicity jets. However, the cuts have been tightened to 
improve the background rejection at higher energies. The efficiency of the selection cuts is 
typically ~33% for inclusive events and ~48% for non-radiative events. The corresponding 
purities of the selected samples are ~88% and ~92% respectively. 

Tracks and electromagnetic clusters, each treated as separate particles with no attempt to 
correct for double-counting of energy, were combined into jets in the following way. First the 
highest energy particle in the event was selected and a cone with a half angle of 35° was defined 
around it. The particle with the next highest energy inside the cone was combined with the 
first. The momenta of the combined particles were added and the direction of the sum was 
used to define a new cone, inside which the next highest energy particle was again sought. This 
procedure was repeated until no more particles were found inside the cone. Similarly, starting 
with the highest energy particle among the remainder, a new cone was initiated and treated 
in the same way. This process continued until finally all the particles in the event had been 
assigned to a cone. 

The following cuts were applied to select r+r" candidates. 

• Hadronic events were rejected by demanding low multiplicity: the number of tracks re- 
constructed in the central tracking detectors had to be at least two and at most six, and 
the sum of the number of tracks plus the number of electromagnetic clusters not more 
than 15. 

• The total energy of an event was restricted in order to reject events from e"'"e~ e"'"e~(7), 
A*^A*~(7) s-iid two-photon processes. The total event energy, E'tot, defined as the scalar 
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sum of all track momenta plus all electromagnetic calorimeter energy, was required to be 

less than l.l^/s. The total electromagnetic calorimeter energy was required to be between 
0.02y^ and 0.7^/s and the scalar sum of track momenta less than 0.8y^. Either the total 
electromagnetic calorimeter energy or the scalar sum of track momenta was required to be 
greater than 0.2^/s. The distribution of Etot/y/s, after all other cuts have been applied, 
is shown in Fig. 7(a) for all centre-of-mass energies combined. The agreement between 
data and simulation is good in the region dominated by the r+r" signal, but poor in 
regions dominated by background; this discrepancy is used to estimate the systematic 
uncertainty in the background. 

Background from two-photon events was further reduced by cuts on the missing momen- 
tum and its direction. The missing momentum in the plane transverse to the beam axis, 
calculated using the electromagnetic calorimeter, was required to exceed 0.015i/i, and 
the polar angle of the missing momentum was required to satisfy | cos^| < 0.99. Fig. 7(b) 
shows the distribution of the missing momentum after all other cuts have been applied, 
for all centre-of-mass energies combined. 

Vertex and TOF cuts were imposed to remove cosmic ray events, as for events. In 

addition, events identified as e+e~ — > ii'^/i" using the criteria described in Section 4.3 
were removed. 

Cones formed from tracks and clusters as described above were classified as either charged 
or neutral. A charged cone was required to contain at least one charged particle with 
transverse momentum greater than 100 MeV and one electromagnetic cluster with energy 
greater than 100 MeV, and the sum of the energy in the electromagnetic calorimeter and 
the track momenta in the cone had to be more than 1% of the beam energy. Neutral 
cones were required to contain no charged particle and an energy in the electromagnetic 
calorimeter of at least 1% of the beam energy. Cones failing these criteria were discarded. 
Events which had exactly two charged cones were selected as e+e~ t~^t~ candidates. 
The direction of each r was approximated by that of the total momentum vector of its 
cone of particles. Events were accepted if both cones satisfied | cos^| < 0.9. To remove 
events with poor momentum reconstruction, the event was rejected if the azimuthal angle 
of either cone, determined using tracks only, lay within 0.5° of an anode plane of the 
central tracking chamber. In addition to the two charged cones, an event may contain 
any number of neutral cones. 

To suppress electron- and muon-pair events further, we reject events with cone energies 
or momenta compatible with these final states. Assuming that the final state consists 
only of two leptons plus a single unobserved photon along the beam direction, the values 
of the polar angles of the two r cones were used to calculate the expected energy of each 
lepton Xi,X2. It was required that 

0.02 < ^{Ef + El)/{X^ + Xl) < 0.8, 

and 

^{P? + Pi)/{X! + Xl) < 0.8, 

where Ei, E2 and Pi, P2 are the total electromagnetic calorimeter energies and scalar sums 
of track momenta, respectively, in each r cone. 
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• Remaining background from e+e" — > e''"e~(7) and e"'"e~e"'"e~ events was reduced by re- 
jecting events if the ratio of the electromagnetic energy to the track momentum in each 
of the r cones was consistent with that expected for an electron. 

• Most of the remaining background from two-photon processes was rejected by a cut on 
the acoUinearity and acoplanarity angles of the two r cones: the acoUinearity angle, in 
degrees, was required to satisfy 

e^oi < (180° - 2tan-^(2mzVi/(s - m|))) + 10° 

and the acoplanarity angle was required to be less than 30°. The value of the cut on 
acoUinearity was chosen such as to include the peak from radiative return events at each 
energy; it is 92° at 200 GeV. The acoplanarity cut was not apphed to events with a photon 
observed in the detector with energy above 30 GeV, if the event was planar (i.e. the sum 
of the opening angles between the three particles was greater than 358°). 

• Events classified as W+W" candidates according to the criteria in [33] were rejected. 

• After the above cuts, the region of y/s' between the radiative return and full energy peaks, 
110 GeV < y/s' < 0.85y^, still contains a significant fraction of background from two- 
photon events. To reduce this background, a hkelihood for the process e+e~ — > t'^t" 
was formed from four variables: the missing momentum calculated using electromagnetic 
clusters, the scalar sum of the track momenta, the invariant mass of the two r cones and 
the difference between the electromagnetic calorimeter energy in the two r cones. The 
value of this likelihood was required to be greater than 0.5. 

• The effective centre-of-mass energy was determined in an identical manner to the 
determination for muon pairs. The observed distribution of y/s' for all energies combined 
is shown in Fig. 9(c). Inclusive events were required to satisfy s'/s > 0.01 and non- 
radiative events were required to have s'/s > 0.7225. 

The selection efficiencies and feedthrough of events from lower s' into the non-radiative sam- 
ples were determined from Monte Carlo events generated with KK2f without interference be- 
tween initial- and final-state radiation, corrected for the four-fermion contribution as discussed 
in Section 2.2. Backgrounds were also determined from Monte Carlo simulations. Efficiencies 
and backgrounds at each energy are summarized in Table 4. 

For the measurement of the angular distributions and asymmetries, the small (~2%) fraction 
of events where the two r cones have the same charge (as determined from the sum of the 
charges of the tracks in the cone) was not used. The final values were obtained by averaging the 
distribution measured using the negative r with that using the positive r, as for the muon pairs. 
The forward-backward asymmetries at each energy were obtained by counting the numbers of 
events in the forward and backward hemispheres, after correcting for background and efficiency. 
The asymmetries were corrected to the full angular range by applying a multiplicative correction 
obtained from Zfitter to the asymmetry measured within the acceptance of the selection cuts 
(Icos^l < 0.9). 
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4.4.2 Systematic Uncertainties 



Systematic errors on the t~^t cross-sections are summarized in Table 11, with a detailed 
breakdown at 200 GeV given in Table 14. The main contributions are discussed below. 

EfRciency. The systematic error on the efficiency was evaluated using high statistics LEP 1 
data and Monte Carlo samples, as for the muon pairs. The t+t~ cross-section at the Z peak 
has been measured with a systematic uncertainty of about 0.5% [39]. As in the case of muon 
pairs, a statistically significant difference between the observed and expected numbers of events 
was seen, and the difference was assigned as the systematic error associated with the selection 
cuts. 

Initial-state radiation modelling. The systematic error on efficiency derived from the 
LEP 1 data does not include the effect of uncertainties in the modelling of initial-state radiation. 
As for hadronic events and muon pairs, this was estimated by reweighting KK2f events from 
0{a'^) to 0{a) CEEX and taking half the predicted change in efficiency. 

Feedthrough. The uncertainty in the feedthrough of events with lower s' into the s' /s > 
0.7225 sample was estimated by comparing the prediction of KK2f with that of KoralZ. 

Backgrounds. The largest background in the tau pairs arises from Bhabha events. Other 
important backgrounds arise from e+e~e"'"e~ and e+e^r+r^ final states. As for the muon pairs, 
systematic errors on each background channel were assessed by comparing data and Monte 
Carlo distributions of selection variables, after loosening selection cuts, in a region enriched in 
the particular background under study. For example, the e'^e" /i'^ /i' background was studied 
using the distribution of total event energy. For small backgrounds which cannot be studied in 
this way, we conservatively assume an error of 50%. 

Interference. The uncertainty arising from the removal of the contribution from interfer- 
ence between initial- and final-state radiation was estimated as described in Section 2.1. 

Asymmetry. Systematic errors on the asymmetry measurement were assessed by compar- 
ing different methods of determining the asymmetry: using tracks, electromagnetic clusters or 
both to determine the r angles. 

4.5 Electron Pairs 

The production of electron pairs is dominated by t-channcl photon exchange, for which a 
definition of s' as for the other channels is less meaningful. In addition, the increased probability 
for final-state radiation relative to initial-state radiation renders the separation between initial- 
and final-state photons more difficult. Events with little radiation were therefore selected by a 
cut on ^acob the acoUinearity angle between electron and positron. We measure cross-sections 
for three different acceptance regions, defined in terms of the angle of the electron, Og-, or 
positron, 9^+, with respect to the incoming electron direction, and the acoUinearity angle: 

• A: I cos^^c-| < 0.9, I cos^e+l < 0.9, ^acoi < 170°; this is a loose 'inclusive' measurement; 

• B: I cos^c- I < 0.7, ^acoi < 10°; this acceptance region is enriched in the s-channel contri- 
bution, and is used for asymmetry measurements; 

• C: I cos^c-| < 0.96, | cos 6*2+ 1 < 0.96, ^acoi < 10°; this 'large acceptance' region is enriched 
in the t-channel contribution and acts as a check on the luminosity measurements. 
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In addition, we measure the electron angular distribution in the region: 

• D: I cos ^e- I < 0.9, ^acol < 10°. 

In all cases, measurements are corrected to correspond to electron and positron energies each 
greater than 0.2 GeV. 

4.5.1 Event Selection 

The selection of e"'"e~ events is identical to previous analyses [1-3]. The selection efficiencies 
are typically ~98%, and the purities of the selected samples ~98%. 

• Events were required to have at least two and not more than eight clusters in the electro- 
magnetic calorimeter, and not more than eight tracks in the central tracking chambers. 

• At least two clusters were required to have an energy exceeding 20% of the beam energy, 
and the total energy deposited in the electromagnetic calorimeter was required to be 
at least 50% of the centre-of-mass energy. For the large acceptance selection, C, which 
has no requirement on the association of tracks to clusters, the total electromagnetic 
energy was required to be at least 70% of the centre-of-mass energy. Distributions of 
total electromagnetic calorimeter energy, after all other cuts, are shown in Fig. 8(b) 
and (c) for acceptance regions B and C for the data from all years combined. There 
is reasonable agreement between data and Monte Carlo. The apparent shght excess of 
data over Monte Carlo at about 80% of the centre-of-mass energy in acceptance B results 
from poor modelling of the energy resolution in the region 0.6 < |cos^| < 0.7. The 
degraded energy resolution in acceptance region C arises from the increased amount of 
material in front of the electromagnetic calorimeter at large | cos^|, where the events are 
concentrated. The detailed modeUing of the electromagnetic calorimeter energy resolution 
has very httle effect on the selection efficiency for e+e~ events. 

• For selections A, B and D, at least two of the three highest energy clusters were required to 
have an associated central detector track. If a cluster had more than one associated track, 
the one with the highest momentum was chosen. If all three clusters had an associated 
track, the two highest energy clusters were chosen to be the electron and positron. For 
the large acceptance selection, C, no requirement was placed on the association of tracks 
to clusters. 

• For the measurement of the forward-backward asymmetry and the angular distribution, 
the two tracks were required to have opposite charge. This extra requirement reduces the 
efficiency by about 3.5% in the region | cos 6'| < 0.9. In addition, due to the extreme charge 
asymmetry for electrons in the forward direction, the problem of charge misassignment 
becomes severe for backward events at small angles. In the measurement of the angular 
distribution we therefore demanded that events with cos^^c- < —0.8 satisfy two extra 
criteria: both electron and positron tracks must have momentum greater than 25% of the 
beam momentum, and there must be only one good track associated with each cluster. 
These criteria significantly reduce the problem of charge misassignment, reducing the 
contamination from wrong-sign events in this region from around 35% to about 15%. 
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• Acceptance cuts on acoUinearity and cos^ were made using the calorimeter clusters, 

with angles corrected for the position of the primary vertex. The acoUinearity angle 
distribution for the inclusive selection, A, is shown in Fig. 8(a), and we see good agreement 
between data and Monte Carlo expectation, including the peak corresponding to s-channel 
radiative return to the Z. 

These cuts have a very high efficiency for e"'"e~ events while providing excellent rejection of 
backgrounds, which either have high multiplicity or lower energy deposited in the electromag- 
netic calorimeter. The efficiency of the selection cuts, and small acceptance corrections, have 
been determined using Monte Carlo events generated with the Bhwide [21] program. These 
were found to be almost independent of energy over the range considered here. Small correc- 
tions have been applied to the efficiencies derived from Monte Carlo simulations to account 
for tracking losses near the central jet chamber anode planes (~ 0.8%), and, in the case of 
the angular distribution, to account for a discrepancy between data and Monte Carlo in the 
fraction of events where both tracks have the same charge (~ 0.5%). Remaining backgrounds 
arise from t^t" events and, in the case of the loose acoUinearity cut, also from electron pairs 
in two-photon events and from radiative Bhabha scattering events in which one electron is 
outside the detector acceptance but the photon is within the acceptance. In the case of the 
large acceptance selection, C, which does not require tracks, the main background arises from 
77 final states. The efficiencies and backgrounds at each energy are summarized in Table 4. 

The forward-backward asymmetries for the 6'acoi < 10° sample at each energy within the 
angular range | cos ^e- 1 < 0.7 were evaluated by counting the numbers of events in the forward 
and backward cos 9^- hemispheres, after correcting for background and efficiency. For both the 
asymmetry and angular distribution measurements, the positive or negative track was used on 
alternate events to reduce systematic effects. 

In Fig. 9(d) we show the distribution of for the inclusive c+e" events for all energies 
combined. The value of s' for each event was estimated from the polar angles of the two electrons 
assuming massless three-body kinematics to calculate the energy of a possible undetected initial- 
state photon along the beam direction as shown in Equation (4). For e'^e", s' is not really well- 
defined, but this calculation gives an estimate of s' for that part of the cross-section proceeding 
via the s-channel. Due to the dominance of the t-channel contributions, for electrons, in contrast 
to the other final states, the radiative return peak forms only a very small contribution. 

4.5.2 Systematic Uncertainties 

The systematic errors associated with the e"'"e~ measurements have generally been estimated 
in a similar manner to the previous analysis [1]. They are summarized in Table 11, with 
a detailed breakdown at 200 GeV given in Table 15. The most significant change is that 
the systematic errors on the differential cross-section measurements (acceptance D) have been 
estimated separately for three cos regions, namely cos 6^- < —0.7, | cos 0^- \ < 0.7 and cos 9^- > 
+0.7, rather than considering the whole distribution together. The systematic errors are not 
expected to be strongly dependent on centre-of-mass energy, so in general they have not been 
estimated separately for each energy point. The most important ones are discussed below. 

Four-fermion contribution. The full size of the change in efficiency arising from includ- 
ing s-channel four-fermion events in the signal definition was included as a systematic error. 
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This affects tfie inclusive selection (selection A) only, and is negligible for events with a tight 
acoUinearity cut. 

Multiplicity cuts. The uncertainties arising from the requirement of low multiplicity 
have been estimated from the change in the number of selected events in data when varying 
the multiplicity cuts used by ±1 unit. 

Calorimeter energy scale and resolution. A detailed comparison between data and 
Monte Carlo has been made of the energy scale and resolution of the electromagnetic calorime- 
ter, and the results of this study used to assess possible effects on the selection efficiency. 
Typically the energy scale was varied by 0.3% and the resolution by 10% of its value. 

Track requirements. Matching between tracks and clusters has been studied using events 
passing all selection cuts, except that only one of the three highest energy clusters has an 
associated track. These events are expected to be mainly e"'"e~7 final states where one electron 
and the photon lie within the acceptance and 77 final states where one photon has converted in 
the detector, with small contributions from other final states. An excess of such events was seen 
in data compared with Monte Carlo expectation. Part of this excess is concentrated in regions 
of (f) near the anode planes of the central jet chamber, and arises from track reconstruction 
problems in this region. The rest could arise from track reconstruction problems, or could 
arise from problems modelling e"'"e~7 or 77 events. For each acceptance region we take the 
excess seen around the jet chamber anode planes plus half the difference between data and 
Monte Carlo in the remaining region of as a correction to the efficiency. This correction is 
typically around 0.8%. Half the difference between data and Monte Carlo in the regions of 
(f) away from the jet chamber anode planes is taken as the systematic error associated with 
track reconstruction. In the case of the angular distribution, the two tracks in an event are 
required to have opposite charge. The fraction of same-sign events in data is roughly 0.5% 
greater than in Monte Carlo. This difference is apphed as a correction to the efficiency derived 
from Monte Carlo, and the value of the correction is included as a systematic error. The higher 
probability of same-sign events in the data is also used to calculate a correction of (1.6±0.5)% 
to the angular distribution in the region cos^^c- < —0.7 arising from charge misassignment. 

Acceptance correction. Because of the steepness of the angular distribution, uncertain- 
ties in the determination of 9 are an important systematic error. These have been assessed 
by comparing measurements of 9 in the electromagnetic calorimeter with those in the central 
tracking chambers and the muon chambers, using e"'"e~ or /x^/i" events as appropriate. These 
studies indicate a possible bias in the 9 reconstruction of electromagnetic clusters of ~1 mrad 
in the endcap region of the detector. The effect of the observed biases on the acceptance was 
calculated using Monte Carlo events, and assigned as a systematic error associated with the 
acceptance correction. 

Background. If a tight acoUinearity cut is applied, the dominant background in the 
selections including tracks is from t'^t~ events. With a loose acoUinearity cut, e''"e~7 and 
e"'"e~e+e~ events are also significant. The systematic error arising from uncertainty in the 
background has been assessed by comparing the numbers of events in data and Monte Carlo 
which pass all cuts except the cut on total calorimeter energy; these events are predominantly 
background. In each acceptance region the larger of the difference between data and Monte 
Carlo or the statistical precision of the test was taken as the associated systematic error. For 
the selection which does not use tracks, acceptance C, the only important background is from 
77 final states; here we used the statistical precision of the OPAL e"'"e~ — > 77 cross-section 
measurement [40] to estimate the uncertainty in this background. 
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Asymmetry. Systematic uncertainties in the asymmetry measurement arise from the 
effects of 9 mismeasurement, charge misassignment and background and efficiency corrections, 
and amount to 0.004. 

4.6 Results 

The numbers of selected events and measured cross-sections for all channels are summarized 
in Table 5. Asymmetries for the leptonic final states are summarized in Table 6, while the 
measured differential cross-sections are given in Tables 7, 8, 9 and 10 for hadrons, r+r" 
and e''"e~ respectively. 

5 Comparison with Standard Model Predictions 

The cross-section and asymmetry measurements at 189-207 GeV are compared with the Stan- 
dard Model predictions in Tables 5-6. Figures 10-13 show cross-sections, for both inclusive 
and non-radiative events, as a function of y^, while Fig. 14 shows the measured asymmetry 
values. The Standard Model predictions are calculated using Bhwide [21] for the e"'"e~ final 
state and Zfitter [6] for all other final states; in this paper we use Zfitter version 6.30 with 
the following input parameters: mz = 91.1852 GeV [39], mtop = 174.3 GeV [41], mniggs = 
115 GeV, Aa^Ji = 0.02761 [42] and Q;s(m|) = 0.1185 [41]. The theoretical uncertainties on the 
cross-section predictions are estimated to be 0.26% for hadronic final states, 0.4% for muon and 
tau final states, 0.5% for electron final states in the endcap region and 2.0% for electron final 
states in the barrel region [7]. In the fits described in Section 7 we assign these values as the 
theoretical errors on the Standard Model cross-sections. For the non-radiative asymmetry val- 
ues we use a theoretical error of 0.004, derived from comparison of the predictions of Zfitter 
and KK2f. The agreement between the measured cross-sections and asymmetry values and the 
Standard Model predictions is generally good. 

The measured differential cross-sections at each energy are given in Tables 7-10. The 
luminosity-weighted averages of data at all energies are compared with Standard Model pre- 
dictions in Figs. 15 and 16. The data arc well-described by the Standard Model curves. 

In order to make a more quantitative test of the compatibility of our cross-sections and 
asymmetries with the Standard Model, we calculate a value between the measurements 
and the Standard Model predictions taking into account statistical and systematic errors and 
their correlations. Correlations between hadron and lepton cross-sections are very small, aris- 
ing mainly from the common luminosity measurements. Correlations between cross-sections 
at different energies for the same channel arise from the systematic uncertainties in both effi- 
ciency and background, but amount at most to 7.6% for hadrons, 4.8% for ii'^ ii~ and 11.4% 
for T'^T~ . Correlations between cross-section and asymmetry measurements are generally neg- 
ligible, amounting at most to about 1% for r'^r", arising from uncertainties in the (mainly 
Bhabha) background. 

The values for the hadronic cross-sections, and the ix^ ix~ and t^t~ cross-sections and 
asymmetries, are shown in Table 16. Note that the 'non-radiative' samples with s'/s > 0.7225 
are a subset of the inclusive events with s' / s > 0.01, so the two x^ values are not independent. 
All measurements are in agreement with the Standard Model expectations. 
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The test would not necessarily reveal a discrepancy in the overall scale of the cross- 
sections or asymmetries compared with the Standard Model expectations. Therefore, as a 
further check, we have also calculated the average value of the ratio of the measurement to 
the Standard Model prediction, using a minimization technique including the experimental 
systematic errors. The results are shown in Table 16. All mean values are compatible with 
unity. The data are thus shown to be compatible with the Standard Model expectations to a 
precision of 1% for hadrons and ~3% for leptons. 

The cross-sections for e"'"e~ events are dominated by the large cos 9 region. Rather than 
comparing the measured integrated cross-sections with the Standard Model, we have calculated 
a for the differential cross-sections, as presented in Table 10. Correlations between cos 9 
bins and between energies are less than 10%, except for the region cos 6* > 0.7 where the 
systematic error is a significant fraction of the total error; in this region the correlation between 
measurements at different energies is 30%-40%. We find a value of 83.3 for 105 degrees of 
freedom, showing excellent agreement. 

5.1 Energy Dependence of agm 

In [1-3] we used non-radiative cross-section and asymmetry measurements to determine the 
electromagnetic coupling constant ttem at LEP 2 energies. We have repeated this fit including 
the new measurements of hadronic, ii^ ijT and t^t~ cross-sections and the [x^ [x~ and t^t~ 
asymmetry values for s' /s > 0.7225 presented here. As before, we form the between the 
measured values and the Standard Model predictions calculated as a function of acm(\/s) using 
Zfitter, with all other Zfitter input parameters fixed [6]. Correlations between measure- 
ments are fully taken into account. In Table 17 we show the results of these fits. We perform 
fits to the data at each energy and also perform a fit to data at all centre-of-mass energies in 
which CKem Tuns with energy with a slope^ obtained from fixing 1/Q;em(0) = 137.036. As input to 
the combined fit we use the new measurements presented here together with the corresponding 
measurements at 130-183 GeV from [2,3]. For the combined fit the value of a^m is quoted at 
the centre-of-mass energy corresponding to the luminosity- weighted average of 1/ s. The errors 
on the fitted values of ccem arise from the errors on the measurements; errors due to uncertain- 
ties in the Zfitter input parameters are negligible. The measured values of are shown in 
Fig. 17. They are consistent with the Standard Model expectations. 

The fits described above use measurements of cross-sections which depend on the measure- 
ment of luminosity. The luminosity measurement assumes the Standard Model running of a^m 
from = to typically = (4 GeV)^, where^ l/ftcm — 134. The fits therefore measure the 
running of ctem only from Qiumi — 4 GeV upwards. To become independent of the luminosity 
measurement, we have repeated the fits replacing the cross-sections for hadrons, muon and tau 
pairs with the ratios a{iiii) / a{qq) and (T(Tr)/(7(qq). This is possible since, above the Z peak, 
hadrons and leptons have very different sensitivity to as discussed in [3]. The results of 
these fits are also shown in Table 17. The values of l/acm are close to those obtained from 
the cross-section fits but with somewhat larger errors. The value of l/ctem obtained from the 
combined fit is 1/Q;em(193.2 GeV) = 126.7l2;|. This is about 4.3 standard deviations below the 
low energy hmit of 137.03599976(50) [45], thus demonstrating the running of ccem from Q"^ — 

^If the Standard Model running of aom is given by — Q;cm(0)/(1 — I1{Q)), then we determine a constant 
K close to 1 such that 0;^*^ = acm(0)/(l — kII(Q}), i.e. the slope dacm/dlnQ is multiplied by k. 

^The hadronic vacuum polarization contribution used in the luminosity measurement is from [46]. 
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to LEP2 energies. This measurement of does not depend on calculations of low-mass 

hadronic loops and is nearly independent of the mass of the Higgs boson and as] it can be 
scaled to the mass of the Z, giving 1/Q;em(91.19 GeV) = 127.9l|8, good agreement with the 
Standard Model prediction of 128.936±0.046 [42]. 

6 S-matrix Analysis 

6.1 Introduction 

Fermion-pair production cross-sections and asymmetries at LEP 1 provide precise informa- 
tion about the Z resonance. The resonance can be described with five parameters: the Z 
mass mz, the Z width Fz, the total hadronic cross-section at the peak, its ratio to the lep- 
tonic cross-section and the leptonic forward-backward asymmetry at the peak. Fitting the 
OPAL data for just these five parameters leads to a precise determination of the Z mass 
mz = 91.1852 ± 0.0030 GeV [39]. However, this fit assumes the contribution to the hadronic 
cross-section from 7 — Z interference behaves as predicted by the Standard Model. A more 
model-independent description of the Z lineshape is provided by the S-matrix approach dis- 
cussed in this section. 

The S-matrix formalism [5] describes the process e+e" — > Z/7* — > ff assuming only the ex- 
change of a combination of two neutral spin-1 bosons of which one is masslcss. Contributions 
from boson exchange and interference are explicitly allowed to vary independently. The result- 
ing parameterizations of the fermion-pair cross-section <J^^^^{s) and asymmetry A^{s) in lowest 
order are: 



Besides the Z mass mz and width Fz there are six parameters per final state fermion f , three for 
the cross-sections and three for the asymmetries. The photon exchange is described by g^"^ and 
gf^{= 0) and is assumed to be known. This leaves four parameters, namely r^"', rp describing 
the Z exchange and jf' for the interference. For the hadronic final state, the parameters 
are summed over all colours and open flavours. Since the hadronic asymmetry is not measured, 
and j'^^j cannot be determined. The lowest order expressions in Equation (5) serve to 
introduce the S-matrix parameters, but cannot be used directly to fit the data without the 
inclusion of large QED radiative corrections. To fit the data, expectations for cross-sections 
and asymmetries depending on S-matrix parameters including QED radiative corrections were 
calculated using the program Smatas Y [47] together with Zfitter [6] . These calculations also 
include very small electroweak corrections to the photon couplings. 

In Equation (5) the Z resonance is described with an s-independent width. UsuaUy loop 
corrections to the Z propagator are absorbed in an s-dependent width via the transformation 
Fz — > sFz/m|. This results in a redefinition of the Z mass and width leading to a numerical 
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shift of: 



mz = mz/\/l + r|/m| ^ mz - 34 MeV , 

Tz = Fz/^l + r|/m| « Tz - 1 MeV . (6) 

Although the S-Matrix parameters are friz and Fz, in this paper the above relations are used 
to give numerical results for mz and Fz to facilitate comparisons with other measurements. 

A fit to LEP 1 data alone leaving the parameters describing the interference free leads to a 
large uncertainty on the Z mass. This is because most LEP 1 data (about 88%) are hadronic 
events taken at three energy points. This results in effectively three very precise measurements 
which dominate the determination of the Z properties. From these three measurements the 
four parameters mz, Fz, r{;°*(j and cannot be determined simultaneously. The interference 
increases the cross-sections at energies above the peak and decreases them at lower energies. 
Therefore a change in this contribution effectively shifts the position of the peak, which can be 
interpreted as a change in the Z mass. This leads to a strong anti-correlation of 96% between the 
fit resuhs of mz = 91.1901 ± 0.0115 GcV and = 0.010 ± 0.650 from LEP 1 data alone [39]. 

LEP 2 data provide additional independent measurements to constrain the contribution 
from interference. Using all OPAL measurements therefore leads to a determination of the Z 
mass which is less model-dependent than that from the five parameter fit to LEP 1 data, with 
an error that is only slightly larger. A similar analysis has been performed previously by other 
experiments [48,49]. Results presented here supersede the OPAL analysis at lower energies [3]. 

6.2 Fit Results 

To derive results for the S-matrix parameters a fit is performed comparing the predictions with 
OPAL measurements of fermion-pair production cross-sections and asymmetries at all LEP 
energies. 

The LEP 1 measurements are described in [39]. The results for the hadronic, e"'"e~, /x^/x" 
and r+r" cross-sections are given in Tables 8-12 of that paper. Leptonic asymmetries can be 
found in Tables 22-24. The measurements used in this analysis are those already corrected for 
the beam energy spread. In the analytical program used to calculate the S-matrix predictions 
t-channel exchange is not implemented. To use Bhabha cross-sections and asymmetries the 
t-channel contribution is corrected in the same way as described in [39]. Alibaba [50] is used 
to determine the Bhabha forward and backward cross-section for the full s + t-channel and for 
the s-channel only. During the fit the difference is added to the S-matrix predictions. The 
t-channel correction is parameterized as a function of (y/s — mz) and thus depends on the value 
of the fitted Z mass. The systematic errors and their correlations are taken as described in [39]. 
The treatment of the errors on the centre-of-mass energy and beam-energy spread is however 
simplified. For the fits described in [39] the effects of these errors were determined iteratively 
during the fitting procedure. In this analysis the effect is calculated prior to the fit according 
to the lineshape determined in [39]. With this simplified error treatment the S-matrix result 
given in [39] can be reproduced, with deviations of all fit parameters less than 3% of their total 
error. 

The high energy measurements used are the non-radiative cross-sections and asymmetries 
for hadrons, /x^/x" and r"^r^ presented here, together with the corresponding results from 
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130-183 GeV [2,3]^. Bhabha cross-sections and asymmetries are not included since even with 
a tight selection^ the observed total cross-section, including s + t-channel, is about an order 
of magnitude larger than the s-channel contribution alone. Correlations between the LEP2 
measurements are taken into account as discussed in Section 5. The systematic errors on the 
LEP 1 and LEP 2 data are basically uncorrelated. The only common error is the systematic 
and theoretical error on the luminosity determination. The resulting effect is found to be small 
and has been neglected. 

Expectations for cross-sections and asymmetries depending on S-matrix parameters are 
calculated using the program Smatasy [47] together with Zfitter [6]. A is calculated 
between the predictions and the measurements of cross-sections and asymmetries. The results 
of fits to OPAL data taken at all LEP energies, with and without the assumption of lepton 
universality, are given in Table 18. The correlation matrices are given in Tables 19 and 20. In 
principle the results for the S-matrix parameters do not depend on the values of the Standard 
Model parameters. However, there is a small effect due to the correction from the s-independent 
width parameters TrTz, Tz to the usual parameters niz, Pz since the width Pz depends on the 
Standard Model parameters. This effect is negligible compared with the overall errors on the 
fitted parameters. As a cross-check a fit is performed using only the LEP 2 data with the 
constraint of the S-matrix result given in [39]. This leads to results very similar to the full fit, 
with differences of all fit parameters being less than 2% of their error. 

As can be seen in Table 18 the fitted Z mass gets smaller by 1.0 MeV when lepton universality 
is imposed. This shift was already observed in a C-parameter fit described in Section 11.1 of [39], 
although it is larger in the S-matrix fit presented here. It is due to a subtle effect in the Bhabha 
i-channel correction, which gives an additional weak constraint on mz (see [39] Appendix B). 

Figure 18 shows the correlation between mz and expressed as confidence level contours 
from fits with lepton universality imposed. Using only LEP 1 data the error on is large 
with a correlation coefficient of -0.96 to mz- Including the LEP 2 data reduces the error on 
ihad ^ factor of five and the correlation is reduced to -0.39. As a result the error on mz is 
much improved from 11.5 MeV (LEP 1 only) to 3.3 MeV. This is now comparable to the error 
of 3.0 MeV obtained from the five parameter fit at LEP 1 which assumes the 7 — Z interference 
according to the Standard Model. Since the fit result for is close to the Standard Model 
expectation the central values for mz are in good agreement. 

7 Limits on New Physics 

New physics could be revealed by deviations of the measured data from Standard Model predic- 
tions. The generally good agreement seen between data and the Standard Model places severe 
constraints on the energy scale of such new phenomena. The new data presented here have 
been combined with previous measurements in order to provide updated limits on four-fermion 
contact interactions. In addition we present limits on the mass of a possible Z' boson. 

■''At ccntrc-of-mass energies of 161 GeV and 172 GeV the cut is s'/s > 0.8. 
^AcoUinearity ^acoi < 10° and angular range for the electron | cos6'e- | < 0.7. 
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7.1 Limits on Four-fermion Contact Interactions 

In the context of composite models of leptons and quarks, a four-fermion contact interaction 
arises as a remnant of the binding force between the substructure of fermions. Alternatively, a 
four-fermion contact interaction could be a good description of deviations from the Standard 
Model due to the exchange of a new very heavy boson of mass mx if irix 3> More generally, 
the contact interaction is considered to be a convenient parameterization to describe possible 
deviations from the Standard Model which may be caused by some unknown new physics. 

In this analysis we consider four-fermion contact interactions which conserve flavour and 
helicity, and in which the SU(3) x SU(2) x U(l) gauge structure of the Standard Model is valid. 
In this framework [51] the Standard Model Lagrangian for e+e~ — > ff is extended by a term 
describing a new effective interaction with an unknown coupling constant g and energy scale 
A: 

where 5 = 1 for e"''e~ — ^ e"''e~ and 6 = otherwise. Here eL(/L) and cr(/r) are chirality projec- 
tions of electron (fermion) spinors, and f]ij describes the chiral structure of the interaction. The 
parameters rjij are free in these models, but typical values are between —1 and depending 
on the type of theory assumed [52]. For example, a coupling of two right-handed currents is 
given by (t^rr = ±1, Vll = ^lr = '7rl = 0). Here we consider the same set of models as in [1-3]. 
The values of rjij which define these models are shown in Table 21. 

The inclusion of a contact interaction modifies both the total cross-section and the angular 
distribution of fermion-pair production. In general, the differential cross-section can be written 
in terms of a parameter e = (gf^/47r)/A^ as 

= ^sm(5, t) + t)s + t)e'' . (8) 

CI COS l7 

Here t = —s{l — cos^)/2 and 6 is the polar angle of the outgoing fermion with respect to the 
e~ beam direction. The C2 term describes the interference between the Standard Model and 
the contact interaction, the C4 term is the pure contact interaction contribution. The exact 
form of these terms depends on the type of fermion in the final state and the particular model 
chosen, and is given, for example, in [53]^. The interference term depends linearly on the rjij 
parameters, and thus can be positive or negative depending on their sign. In fits to the data, 
the Standard Model cross-sections o-sM('5-i) were calculated using Bhwide for the c+c^ final 
state and Zfitter for all other final states. Radiative corrections to the lowest order contact 
interaction terms were taken into account as described in [3]. 

We have fitted the measurements of the non-radiative cross-sections for e"'"e~ — > qq, non- 
radiative cross-sections and asymmetries for e"^e~ — > and e"'"e~ — > t+t~ and the differ- 
ential cross-sections for e+e~ — > e+e~ at 189 GeV to 207 GeV presented here, together with 
the corresponding measurements at 130-183 GeV [2,3]. In all cases we use a fit, including 
the correlated systematic errors between the measurements and theoretical uncertainties in the 
Standard Model predictions as discussed in Section 5. Fits are performed with the parameter 
e as the fitting parameter. The results for positive and negative interference with the Standard 

^Equation (2) in [53] has a typographical error: the factor 4s on the left-hand side should be replaced by 2s. 
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Model (i.e. the sign of the rjij parameters) are equivalent under the transformation e -(-^ —s; it 
is therefore sufficient to fit only for the case of positive interference but to allow e to be both 
positive and negative. Limits on the energy scale A were extracted assuming 5'^/47r = 1. The 
95% confidence limits correspond to a change in of 3.84. 

The results are shown in Table 21 and illustrated graphically in Fig. 19. The limits for 
qq are derived from the hadronic cross-sections assuming the new interaction couples to all 
flavours equally. Those for up-type quarks and down-type quarks are obtained by fitting the 
hadronic cross-sections assuming the new interaction couples only to one fiavour, whereas those 
for uu + dd assume a coupling to one generation only. The combined results include all leptonic 
channels and the hadronic cross-sections. The two sets of values A_|_ and A_ shown in Table 21 
correspond to positive and negative values of e respectively, refiecting the two possible signs 
of rjij in Equation (7). The data are particularly sensitive to the VV and AA models; the 
combined data give lower limits on A in the range 13-16 TeV for these models. For the other 
models the lower limits generally lie in the range 9-13 TeV. The limits are typically 1 TeV 
higher than those for 130-189 GeV data alone [1]. 

Contact interactions involving quarks have also been studied in ep and pp collisions, where 
limits comparable to our values are found [54,55]. Atomic physics parity violation experiments 
can place higher hmits (~ 15 TeV [56]) on models of eeuu and eedd contact interactions which 
violate parity. 



7.2 Limits on a Z' Boson 
7.2.1 Z' Model Predictions 

Many theories predict a second heavy neutral vector boson Z'° in addition to the Standard 
Model gauge boson Z°. The Z° has fermion couplings as predicted by the Standard Model, 
whereas the axial and vector coupling constants of the Z'" to fermions are parameters of the 
particular model. In general the additional heavy boson Z'" will mix with the Z'^ boson. The 
observed particles are the mass eigenstates Z and Z' [57]: 

Z \ _ / cos6'm sin^^M \ 
Z' J ~ [ -sin^M cos^M ; V Z" 

The mixing angle is a free parameter of the model. Here and in the following Z" and Z'" 
denote gauge eigenstates whereas the mass eigenstates formed by mixing within the Z' model 
are denoted as Z and Z'. 

In this paper we consider several Z' models. In E6 GUT, the E(6) group may incorporate 
the Standard Model groups of colour SU(3)c, weak isospin SU(2)l and hypercharge U(1)y in 
the following way [58] : 

E(6) ^ SO(IO) X U(l)^ 
SO(IO) ^ SU(5) X U(l)^ 
SU(5) ^ SU(3)c X SU(2)l x U(1)y • (10) 

Thus the two additional gauge groups U(l);^ and U(l)^ are introduced, each related to a new 
gauge boson. In general, the Z'° will be a mixed state of these two groups: 

Z'° = Z^ sin^E6 + 'Z'x cos6'e6 ■ (11) 
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We derive limits on the Z' mass and mixing angle for all values of 9^6, with particular emphasis 
on three special cases: the two no-mixing models x ^iid ip {6^6 = and 7r/2) and the r) model 
with = — arctany's/S = —0.91. In this latter case the E6 group is broken by a non-Abelian 
discrete symmetry to a rank-5 group which may occur in superstring theories [58] . 

Another approach is the left-right symmetric model (LR) [59]. In this scenario a symme- 
try group SU(2)r is introduced whose three vector bosons couple to right-handed fermions. 
The coupling constants of the Z'° to fermions depend on one parameter cklr. This can take 
values ^2/3 < cklr < ^J (cos^ ~ sin^ ^w) / sin^ where is the weak mixing angle. For 
cclr = 1.53 (the upper limit of the allowed range) left and right-handed coupling constants are 
approximately the same. For q;lr = ^2/3 the LR model is equivalent to the E6 x model. We 
derive limits on the Z' mass and mixing angle as a function of Olr; of particular interest is 
the symmetric case of equal left- and right-handed couplings. The axial and vector coupling 
constants of the Z'° to fermions for the E6 and LR models can be found in [58,60]. 

In addition to the models described above, we also present limits in the case of a sequential 
Standard Model (SSM) Z', which has the same couphngs to fermions as the Standard Model Z. 

In a more model-independent approach, the Z' can be directly described in terms of its axial 
and vector couplings to fermions, a'^ and v'^ [61]. At energies far from the Z' resonance, the 
data are sensitive to the normalized couplings 

N / N / 

at ~ af\l , Vf ^ Vf\ . 

•' •' V 47rmz' V 47rmz' 

We present limits on the leptonic couplings, assuming lepton universality and 5'^/47r = 1. 

In this analysis, cross-sections and asymmetries predicted by the Z' models are obtained 
using Zfitter [6] together with Zefit [60]. Input parameters to the routine are the model 
parameters: masses of Z and Z', model angle 9e6 or q;lr and mixing angle between Z and Z'. 
In addition the usual Zfitter input parameters, as described in Section 5, are used. For the 
Z' mass range relevant here (300 < mz' < 5000 GeV) the difference between the fitted Z mass 
within the Standard Model and the mass within the Z' model is less than 1 MeV, well within 
the experimental uncertainty of the OPAL measurement, mz — 91.1852 ± 0.0030 GeV [39]. In 
the fits described below, mz is treated as a free parameter. The parameters ctg, mtop and mmggs 
are fixed to the values given in Section 5 unless otherwise stated. 

7.2.2 Analysis and Results 

Cross-sections and asymmetries measured at energies around the Z peak give a precise deter- 
mination of the properties of the Z boson. If this particle is not the Standard Model Z" but a 
mixture with Z'^ the couplings to fermions will change. In particular, the measured width of the 
Z is sensitive to the mixing angle ^m, and this angle is therefore constrained by the LEP 1 data. 
At energies above 130 GeV the interference between Z and Z' becomes increasingly important 
and the data are very sensitive to the mass of the Z'. Since changes to cross-sections and 
asymmetries arise from interference terms, the precise form of these changes depends strongly 
on the model. 

To obtain limits on the Z' properties, cross-sections for the processes e+e~ —>■ yU"*"/^", r+r ~, qq 
and the forward-backward asymmetries for the leptonic processes e"'"e~ fi~^fi~ ,t~^t~ at en- 
ergies around the Z resonance [39] and the non-radiative values at ^/s = 130-183 GeV [2, 3] 
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and 189-207 GeV presented here are compared to the predictions of the Z' models. Cross- 
sections and asymmetries for b- and c-quark production arc found to yield negligible additional 
sensitivity and are not used. In calculating the between the predictions of the model and 
the measurements the correlations of the experimental errors are taken into account as in the 
S-matrix fit. 

A between model predictions and the measurements is calculated for different values 
of the Z' mass and the mixing angle ^m- The difference between the minimum and the 
Standard Model is small for all models, being at most 0.8 for the LR model. 

In Fig. 20 we show the 95% confidence level exclusion contours in the rriz' - Om plane for the 

E6 models x^'^-V f-i^*^ the LR model with q;lr = 1.53. These correspond to > Xmm + 5.99. 
The limit on the Z' mass depends strongly on the model. The allowed range for the mixing 
angle is approximately -2 mrad < 6*]^ < 3 mrad, and only for the E6 rj model is the contour 
much broader. For very large Z' masses the contours for the different models become similar. 

In Fig. 21(a) and (b) we present one-dimensional limits on mz' and for a scan over the 
E6 and LR model angles. The 95% confidence level hmits correspond to x^ > Xmm + 3-84. 
The Standard Model parameters were again fixed. Numerical values of the limits on mz' and 

for parameters corresponding to the Xii^i'H ^'^^ symmetric LR models obtained from these 
one-dimensional fits are given in Table 22, where we also present corresponding limits on the 
sequential Standard Model Z'. 

To assess the effect of fixing the Standard Model parameters, we have also performed fits 
in which the strong coupling constant and the top quark mass were treated as free parameters, 
but constrained by their experimental error. The Higgs boson mass was fixed. The sensitivity 
of the LEP 1 data to JTiHiggs arises mainly from the Z width, which changes in Z' models. The 
data cannot discriminate between Z' and Higgs effects, therefore a fit with a free Higgs boson 
mass would be numerically unstable. Figure 21(c) shows the change in the limit on the mixing 
angle as a function of model parameters. The limit on the mixing angle obtained with and 
mtop free is typically 10% less restrictive than that obtained when fixing these parameters. 
Changes to the limits on the Z' mass are at most a few GeV. Since the Higgs boson mass is 
unknown its infiucnce on the limits has been studied by performing fits with a different fixed 
mass. Figure 21(d) shows the change in the limit on the mixing angle when the Higgs boson 
mass is set to 250 GeV rather than its default value of 115 GeV. The change is generally small, 
but amounts to almost 30% in the region of the E6 model. 

Limits on the vector and axial-vector couplings of a Z' boson to leptons within the model- 
independent framework were derived with the mixing angle set to zero. Therefore only 
the leptonic data at energies of 130 GeV and above were used in the fit, and the Z mass 
was fixed. The couplings cannot be determined independently from the Z' mass, so we have 
determined limits for fixed masses of 300, 500 and 1000 GeV. The 95% confidence level limits on 
the vector and axial- vector couplings, v'^ and a[ are shown in Fig. 22. The exclusion contours 
are roughly rectangular in shape. In terms of the normalized couplings, calculated at the 
luminosity-weighted mean centre-of-mass energy of 193.2 GeV, we find 95% confidence level 
limits of |a?| < 0.145 and \vf\ < 0.127. 
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8 Conclusions 



Wc have presented new measurements of cross-sections and asymmetries for hadronic and 
lepton-pair production in e^e~ collisions at centre-of-mass energies between 189 GeV and 
209 GeV. At these energies, strong initial-state radiation leads to excitation of the Z. We 
therefore distinguish two kinematic regions depending on s', the square of the centre-of-mass 
energy of the e+e~ system after initial-state radiation: an 'inclusive' region with s'/s > 0.01 
and a 'non-radiative' region with s'/s > 0.7225. The results for both inclusive fermion-pair 
production and for non-radiative events are in good agreement with Standard Model expecta- 
tions. From these and earlier measurements we derive a value for the electromagnetic coupling 
constant 1/Q;em(193.2 GeV) = 127.4l|o- addition, the results have been used together with 
OPAL measurements at 91-183 GeV within the S-matrix formalism assuming lepton univer- 
sality to determine the 7-Z interference term = 0.144 ± 0.078 and to make an almost 
model-independent measurement of the Z mass, ruz = 91.1872 ± 0.0033 GeV. 

The measurements have also been used to place limits on new physics. In the context of 
a four-fermion contact interaction we have improved the limits on the energy scale A from 
typically 3-13 TeV to 5-16 TeV, assuming g^/4:7i = 1. Lower limits on the mass of a possible 
Z' boson in the range 334 GeV to 1018 GeV, depending on the model, have been obtained. 
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Interference Corrections s' / s > 0.01 


yfs / GeV 


Aa/asM(qq) (%) 


Aa/asuil^l^) (%) 


Aa/asuirr) (%) 


A/lfb(Ai/x) 


AAfb(rr) 


189 


+0.04±0.06 


-0.63±0.01 


-0.47±0.04 


-0.0070±0.0002 


-0.0057±0.0001 


192 


+0.04±0.06 


-0.63±0.01 


-0.49±0.04 


-0.0070±0.0002 


-0.0058±0.0002 


196 


+0.05±0.05 


-0.64±0.01 


-0.47±0.04 


-0.0071±0.0002 


-0.0058±0.0002 


200 


+0.06±0.04 


-0.65±0.01 


-0.48±0.04 


-0.0072±0.0002 


-0.0059±0.0002 


202 


+0.06±0.04 


-0.65±0.01 


-0.51±0.03 


-0.0072±0.0002 


-0.0061±0.0001 


205 


+0.07±0.05 


-0.65±0.01 


-0.51±0.03 


-0.0073±0.0002 


-0.0060±0.0002 


207 


+0.08±0.05 


-0.65±0.01 


-0.50±0.04 


-0.0073±0.0002 


-0.0061i0.0001 


Interference Corrections s' / s > 0.7225 


/ GeV 


A(T/asM(qq) (%) 


Aa/asuil^l^) (%) 


Aa/asuirr) (%) 


AAfbii^n) 


AAfb(rr) 


189 


+0.14±0.24 


-1.46±0.12 


-0.95±0.02 


-0.0133±0.0015 


-0.101±0.0009 


192 


+0.17±0.21 


-1.46±0.12 


-0.96±0.02 


-0.0135±0.0015 


-0.102±0.0009 


196 


+0.22±0.18 


-1.46±0.12 


-0.93±0.02 


-0.0136±0.0015 


-0.101±0.0008 


200 


+0.27±0.13 


-1.46±0.12 


-0.94±0.02 


-0.0137±0.0015 


-0.103±0.0009 


202 


+0.30±0.11 


-1.46±0.12 


-0.96±0.02 


-0.0137±0.0015 


-0.105±0.0009 


205 


+0.34=b0.09 


-1.46±0.12 


-0.96±0.02 


-0.0139±0.0015 


-0.105±0.0009 


207 


+0.37±0.10 


-1.46±0.11 


-0.95±0.02 


-0.0139±0.0015 


-0.105±0.0009 



Table 1: Corrections Act and AAf^ which have been applied to the measured cross-sections and asymmetries in order to remove the 
contribution from interference between initial- and final-state radiation. Cross-section corrections are expressed as a percentage of the 
expected Standard Model cross-sections, calculated using Zfitter, while asymmetry corrections are given as absolute numbers. 



v/i / GeV 






nominal 


range 


mean 


/ Cdt / pb-i 


A£//: (%) 


189 




188.635±0.040 


185 


0.21 


192 




191.590±0.042 


29 


0.32 


196 




195.526it0.042 


77 


0.26 


200 




199.522±0.042 


79 


0.27 


202 




201.636±0.012 


38 


0.30 


205 


202.5 - 205.5 


204.881±0.050 


82 


0.26 


207 


205.5 - 209.0 


206.561±0.050 


137 


0.24 



Table 2: Centre-of-mass energy range, luminosity- weighted mean centre-of-mass energy [16], 
approximate integrated luminosity collected and total error on the luminosity measurement at 
each nominal energy point. The precise amount of data used in each analysis varies slightly 
from channel to channel. 



v/i / GeV 


189 


192 


196 


200 


202 


205 


207 


Experimental systematic 


0.15 


0.20 


0.18 


0.19 


0.19 


0.18 


0.18 


Beam energy 


0.04 


0.04 


0.04 


0.04 


0.04 


0.05 


0.05 


Theory 


0.12 


0.12 


0.12 


0.12 


0.12 


0.12 


0.12 


Data statistics 


0.09 


0.20 


0.13 


0.13 


0.19 


0.13 


0.10 


Total 


0.21 


0.32 


0.26 


0.27 


0.30 


0.26 


0.24 



Table 3: Errors on the luminosity measurement (in %) at each nominal centre-of-mass energy. 
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Efficiencies and backgrounds at ^/s = 189 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq 


s'/s > 0.01 






86.6±0.3 


6.1±0.7 






s'/s > 0.7225 






87.0±0.3 


1.55±0.06 


0.959±0.043 




s'/s > 0.01 






74.7±0.8 


0.42±0.02 






s'/s > 0.7225 






87.9±0.9 


0.058±0.006 


0.044±0.005 




s'/s > 0.01 






33.8±0.8 


0.32±0.03 






s'/s > 0.7225 






48.4±1.2 


0.123±0.014 


0.055±0.003 


c+e^ 


A: |cos6'e±| < 0.9, 


^acol 


< 170° 


97.5±0.4 


1.8±0.3 






B: 1 cos 6*0-1 < 0.7, 


^acol 


< 10° 


99.0±0.3 


0.25±0.03 






C: 1 cotiO^± \ < 0.96, 0ac, 


1 < 10" 


98.5±0.1 


10.i±0.5 




Efficiencies and backgrounds at -y/i = 192 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq 


s'/s > 0.01 






85.9±0.3 


6.2±0.7 






s'/s > 0.7225 






86.7±0.3 


1.54±0.06 


0.972±0.044 


^+//- 


s'/s > 0.01 






74.4±0.8 


0.43±0.02 






s'/s > 0.7225 






87.9±0.9 


0.062±0.005 


0.042±0.004 




s'/s > 0.01 






33.6±0.8 


0.31±0.03 






s'/s > 0.7225 






48.4±1.2 


0.121±0.013 


0.052±0.003 


e+e"" 


A: cos0c± < 0.9, 


^acol 


< 170° 


97.5±0.4 


1.7±0.3 






B: 1 cos6'e-| < 0.7, 


^acol 


< 10° 


98.9±0.3 


0.24±0.02 






C: coseieil < 0.96, 


^acol < 10° 


98.5±0.4 


10.0±0.5 




Efficiencies and backgrounds at -y/i = 196 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq 


s'/s > 0.01 






85.1±0.3 


6.2±0.7 






s'/s > 0.7225 






86.2±0.3 


1.52±0.06 


0.868±0.039 




s'/s > 0.01 






74.0±0.8 


0.45±0.02 






s'/s > 0.7225 






87.9±0.9 


0.068±0.005 


0.040±0.004 




s'/s > 0.01 






33.3±0.8 


0.31±0.03 






s'/s > 0.7225 






48.4±1.1 


0.118±0.012 


0.049±0.003 


e+e" 


A: |cos6'e± < 0.9, 


^acol 


< 170° 


97.4±0.4 


1.7±0.2 






B: 1 cos 0e- 1 < 0.7, 


^acol 


< 10° 


98.9±0.3 


0.23±0.02 






C: 1 cos 6*6*1 < 0.96, 


6'acol < 10° 


98.5±0.4 


9.7±0.5 





Table 4: Efficiency of selection cuts, accepted background, and feedthrough of events generated 
with lower s' into the non-radiative samples, for each channel at each energy. The (very small) 
contribution of events with s'/s < 0.01 to the inclusive sample is included in the efficiency. The 
errors include Monte Carlo statistics and systematic effects. In the case of electron pairs, the 
efficiencies are effective values including the efficiency of selection cuts for events within the 
acceptance region and the effect of acceptance corrections. An acceptance of |cos^c±| < 0.9 
(or 0.96) means that both electron and positron must satisfy this cut, whereas | cos^e-| < 0.7 
means that only the electron need do so. 
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Efficiencies and backgrounds at y/s — 200 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq s' /s > 0.01 

s'/s > 0.7225 


84.2±0.4 
85.8±0.3 


6.2±0.6 
1.50±0.06 


0.832±0.037 


fi+fi- s'/s > 0.01 

s'/s > 0.7225 


73.7±0.8 
87.9±0.9 


0.47±0.02 
0.073±0.005 


0.038±0.004 


r+r- s'/s > 0.01 

s'/s > 0.7225 


32.9±0.8 
48.4±1.1 


0.30±0.02 
0.115±0.011 


0.045±0.002 


e+e~ A 
B 
C 


|cos6»c±| < 0.9, 6'acoi < 170° 

|cOS0e-| <0.7, ^acol < 10° 
|cOs6'e±| < 0.96, 6lacol < 10° 


97.4±0.4 
98.9±0.3 
98.4±0.4 


1.6±0.3 
0.22±0.02 
9.3±0.5 




Efficiencies and backgrounds at ^/s = 202 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


83.8±0.4 

85.6±0.3 


6.2±0.6 

1.49±0.06 


0.819±0.036 


fi^li- s'/s > 0.01 

s'/s > 0.7225 


73.5±0.8 
87.9±0.9 


0.47±0.02 
0.076±0.005 


0.037±0.004 


r+r- s'/s > 0.01 

s'/s > 0.7225 


32.8±0.8 

iS.lii.l 


0.29±0.02 

0.11 i±0. Oil 


0.013±0.002 


e+e~ A 
B 
C 


1 COS 6*3*1 < 0.9, 6'acoi < 170° 

|cOS^e-| < 0.7, 6'acoi < 10° 
|cOs6'e±| < 0.96, 6lacol < 10° 


97.4±0.4 
98.9±0.3 
98.4±0.4 


1.6ib0.3 
0.22±0.02 
9.1±0.5 




Efficiencies and backgrounds at ^/s = 205 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


83.1±0.4 
85.2±0.3 


6.3±0.6 
1.48±0.05 


0.805±0.036 


fi+li~ s'/s > 0.01 

s'/s > 0.7225 


73.2±0.8 
87.9±0.9 


0.49±0.02 
0.080±0.006 


0.035 ±0.004 


r+r" s'/s > 0.01 

s'/s > 0.7225 


32.5±0.8 
48.4±1.1 


0.29±0.02 
O.llliO.OlO 


0.040 ±0.002 


e+e~ A 
B 
C 


|cOs6'e±| < 0.9, 6'acoi < 170° 

|cOS0e-| <0.7, 0acol < 10° 

1 COS 6*6*1 < 0.96, 6'acoi < 10° 


97.2±0.5 

98.8±0.4 
98.4±0.4 


1.5±0.3 

0.21±0.02 
8.8±0.4 




Efficiencies and backgrounds at s/s = 207 GeV 


Channel 


Efficiency (%) 


Background / pb 


Feedthrough / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


82.7±0.4 
85.0±0.3 


6.3±0.6 
1.47±0.05 


0.746±0.033 


fi+ij.- s'/s > 0.01 

s'/s > 0.7225 


73.0±0.8 

87.9±0.9 


0.50±0.03 

0.082±0.006 


0.034 ±0.003 


r+r- s'/s > 0.01 

s'/s > 0.7225 


32.4±0.8 
48.4±1.1 


0.28±0.02 
O.llOzbO.Oll 


0.039 ±0.002 


e+e~ A 
B 
C 


Icoseieil < 0.9, 6'acoi < 170° 
|cOS0e-| <0.7, 0acol<lO° 
1 cos 6'e± 1 < 0.96, 6'acoi < 10° 


97.2±0.5 
98.8±0.4 
98.4±0.4 


1.5±0.3 
0.21±0.02 
8.7±0.4 
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Cross-sections at ^/s = 189 GeV 


Channel 


/£dt/pb-i 


Events 


cr / pb 


cjSM / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


185.9 


17146 
4019 


99.5±0.8±0.9 
22.0±0.4±0.1 


98.9 
22.2 


H+H- s'/s > 0.01 

s'/s > 0.7225 


178.6 


1128 
519 


7.85±0.25±0.09 
3.14±0.15±0.03 


7.75 
3.21 


T+T- s'/s > 0.01 

s'/s > 0.7225 


179.0 


554 
333 


8.17±0.39±0.21 
3.45±0.21±0.09 


7.74 
3.21 


e+e" A 
B 


1 cos 6*0*1 < 0.9, l9acol < 170° 
\cOSe,- \ < 0.7, 6'acol < 10° 

cos^^,,±| < 0.9G, (9„,„i < 10° 


185.9 


20538 
3758 

576(39 


111.5±0.8±0.6 
20.2±0.3±0.1 

301.4±1.3±1.5 


110.3 
20.1 

307.7 


Cross-sections at -y/i = 192 GeV 


Channel 


/£dt/pb"i 


Events 


cr / pb 


cjSM / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


29.6 


2617 
643 


95.9±2.0±0.9 
22.2±0.9±0.1 


95.0 
21.3 


fi+ft- s'/s > 0.01 

s'/s > 0.7225 


29.0 


173 
77 


7.40±0.61±0.09 
2.86±0.34±0.03 


7.47 
3.10 


T+r" s'/s > 0.01 

s'/s > 0.7225 


29.1 


85 
50 


7.74±0.95±0.20 

3.17±0.50±0.08 


7.47 

3.10 


e+e" A 
B 
C 


|cOs6'e±| < 0.9, 6'acol < 170° 
\cOSe,-\ < 0.7, ^acol < 10° 
|cOs6'e±| < 0.96, 6lacol < 10° 


29.5 


3084 
577 
9034 


105.6±1.9±0.6 

19.5±0.8±0.1 

301.0±3.3±1.6 


106.9 
19.5 
298.3 


Cross-sections at -y/s = 196 GeV 


Channel 


jCdt/ph-' 


Events 


a / pb 


(tSM / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


77.8 


6351 
1509 


88.8±1.2±0.9 
19.8±0.6±0.1 


90.2 
20.2 


At+At" s'/s > 0.01 

s'/s > 0.7225 


76.0 


435 
207 


7.08±0.37±0.08 
2.93±0.22±0.03 


7.13 
2.96 


T+r" s'/s > 0.01 

s'/s > 0.7225 


75.9 


206 
120 


7.21±0.57±0.19 
2.89±0.30±0.07 


7.12 
2.96 


e+e" A 
B 
C 


|cOs6'e±| < 0.9, 6'acol < 170° 
\cOSe,-\ < 0.7, ^acol < 10° 
|cOs6'e±| < 0.96, 6'acol < 10° 


77.7 


7879 
1448 
22618 


102.3±1.2±0.5 

18.6±0.5±0.1 

285.7±2.0±1.5 


102.6 
18.7 
286.5 



Table 5: Measured cross-sections, integrated luminosity used in the analysis and numbers of 
selected events at each energy. For the cross-sections, the first error shown is statistical, the 
second systematic. As in [3], the cross-sections for hadrons, /^+/^" and r+r^ are defined to cover 
phase-space up to the limit imposed by the s'/s cut, with y/s' defined as the invariant mass of 
the outgoing two-fermion system before final-state radiation. The contribution of interference 
between initial- and final-state radiation has been removed. The last column shows the Standard 
Model cross-section predictions from Zfitter [6] (hadrons, //+/x", r+r") and Bhwide [21] 
(e+e"). 
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Cross-sections at y/s — 200 GeV 


Channel 


jCdt / 


Events 


cr / pb 


ctSM / 


qq s' /s > 0.01 

s'/s > 0.7225 


79.4 


6100 
1468 


83.9±1.2±0.9 
18.9±0.5±0.1 


85.7 
19.1 


s'/s > 0.01 
s'/s > 0.7225 


78.2 


423 
202 


6.67±0.36±0.08 
2.77±0.21±0.03 


6.80 
2.83 


T+r- s'/s > 0.01 

s'/s > 0.7225 


78.0 


205 
132 


7.04±0.56±0.18 
3.14±0.30±0.08 


6.80 
2.83 


e+e~ A 
B 
C 


|cOs6'e±| < 0.9, 6'acoi < 170° 
\cOSe,-\ < 0.7, 0acol < 10° 
|cOs6'e±| < 0.96, 6'acoi < 10° 


79.4 


7819 
1444 
22046 


99.5±1.1±0.5 
18.2±0.5±0.1 
272.8±1.9±1.4 


98.5 
17.9 
275.2 


Cross-sections at = 202 GeV 


Channel 


jCdt/ph-' 


Events 


cr / pb 


a^^ / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


38.2 


2898 

692 


83.2±1.7±0.9 

18.5±0.8±0.1 


83.5 

18.6 


li+li- s'/s > 0.01 

s'/s > 0.7225 


36.8 


171 
82 


5.63±0.48±0.07 
2.36±0.28±0.03 


6.64 
2.77 


T+r- s'/s > 0.01 

.s-'/.s > 0.7225 


36.9 


104 

59 


7.69±0.84±0.20 

2.95±0.13±0.07 


6.63 
2.77 


e+e" A 
B 
C 


|cos6'e±| < 0.9, 6lacol < 170° 
IcOS^e-l < 0.7, 0acol < 10° 
|cOs6'e±| < 0.96, 6'acoi < 10° 


38.2 


3697 
682 
10551 


97.8±1.6±0.5 
17.8±0.7±0.1 
271.7±2.7±1.4 


96.4 
17.6 
269.5 


Cross-sections at y/s = 205 GeV 


Channel 


JCdt/ph-^ 


Events 


cr / pb 


/ pb 


qq s'/s > 0.01 

s'/s > 0.7225 


82.3 


6094 
1458 


81.7±1.1±0.9 
18.2±0.5±0.1 


80.2 
17.8 


H+H' s'/s > 0.01 

s'/s > 0.7225 


78.9 


418 
212 


6.53±0.35±0.08 
2.88±0.21±0.03 


6.41 
2.67 


r+T" s'/s > 0.01 

s'/s > 0.7225 


78.9 


199 
117 


6.84±0.55±0.18 
2.72±0.28±0.07 


6.40 
2.67 


e+e~ A 
B 
C 


1 COS 6'e± 1 < 0.9, 6'acoi < 170° 
1 COS 0,-1 < 0.7, 0acol < 10° 
|cOs6'e±| < 0.96, 6'acoi < 10° 


82.3 


7613 

1433 
21916 


93.6±1.1±0.6 

17.4±0.5±0.1 
261.8±1.8±1.3 


93.4 
17.0 
261.1 


Cross-sections at ^/s = 207 GeV 


Channel 


/£dt/pb-i 


Events 


cr / pb 


fjSM / pb 


qq s'/s > 0.01 

s'/s > 0.7225 


137.4 


9686 
2260 


77.7±0.9±0.8 
16.8±0.4±0.1 


78.6 
17.5 


s'/s > 0.01 
s'/s > 0.7225 


134.4 


739 

347 


6.81±0.28±0.08 
2.77±0.16±0.03 


6.29 

2.63 


T+r" s'/s > 0.01 

s'/s > 0.7225 


134.5 


318 
203 


6.39±0.41±0.17 
2.78±0.22±0.07 


6.28 
2.63 


e+e- A: | cos e^± \ < 0.9, 6'acoi < 170° 

B: |cOs6'e-| < 0.7, 6'acoi < 10° 

C: cos6'e±| < 0.96, 6'acoi < 10° 


137.7 


12335 
2320 
35997 


90.6±0.8±0.6 
16.9±0.4±0.1 
257.1±1.4±1.3 


91.8 
16.7 
256.9 
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Asymmetries at -y/s 


= 189 GeV 






Nf 


iVb 




/iSM 
^fb 






s'/s > 0.01 


733 


395 


0.252±0.030±0.006 


0.281 






s'/s > 0.7225 


399.5 


119.5 


0.548±0.039±0.005 


0.569 


r^ 




s'/s > 0.01 


373 


170 


0.304±0.046±0.008 


0.281 






s'/s > 0.7225 


253.5 


73.5 


0.591±0.054±0.012 


0.569 






cos6le- < 0.7 

aiKl 0,,,,,i < 10° 


3332 


354 


0.811±0.010±0.004 


0.814 






Asymmetries at 


= 192 GeV 






Nf 


iVb 


^fb 


/iSM 
^fb 




V- 


s'/s > 0.01 


100 


73 


0.095±0.080±0.006 


0.280 






s'/s > 0.7225 


52 


25 


0.341±0.115±0.005 


0.566 


r^ 




s'/s > 0.01 


62 


21 


0.444±0.111±0.008 


0.280 






s'/s > 0.7225 


42.5 


6.5 


0.813±0.109±0.013 


0.565 






cos 6*6- < 0.7 

and ^aroi < 10" 


518 


46 


0.841±0.023±0.004 


0.814 






Asymmetries 


at 


= 196 GeV 






ATf 


A^b 




/iSM 


1^' 


V- 


s'/s > 0.01 


305.5 


129.5 


0.358±0.048±0.005 


0.279 






s'/s > 0.7225 


172.5 


34.5 


0.683±0.055±0.005 


0.562 


r^ 




s'/s > 0.01 


125.5 


72.5 


0.175±0.077±0.008 


0.279 






s'/s > 0.7225 


78 


36 


0.373±0.103±0.013 


0.561 






cos6le- < 0.7 
and e^o\ < 10° 


1276 


137 


0.810±0.016±0.004 


0.815 



Table 6: Measured asymmetry values and numbers of forward (A^f ) and backward (A^b) events at 
each energy. The measured asymmetry values include corrections for background and efficiency, 
and in the case of and t'^t~ are corrected to the full solid angle with interference between 
initial- and final-state radiation subtracted. The first error shown is statistical, the second 
systematic. The final column shows the Standard Model predictions of Bhwide for e'^e" and 
Zfitter for the other final states. 
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Asymmetries at -y/s 


= 200 GeV 












/iSM 
^fb 




s'/s > 0.01 


294.5 


128.5 


0.346±0.049±0.005 


0.278 






s'/s > 0.7225 


164 


38 


0.637±0.059±0.005 


0.558 


T" 




s'/s > 0.01 


142 


59 


0.341±0.074±0.005 


0.279 






s'/s > 0.7225 


107 


23 


0.700±0.077±0.009 


0.558 




'e~ 


1 cos 6^6- 1 < 0.7 
and ^acoi < 10° 


1290 


142 


0.805±0.016±0.004 


0.815 






Asymmetries at y/s 


- 202 GeV 






A^f 


A^b 


-4fb 


/iSM 

ID 






s'/s > 0.01 
s'/s > 0.7225 


114 

61 


57 

21 


0.277±0.080±0.006 

0.189±0.101±0.005 


0.278 

0.556 


r" 




s'/s > 0.01 


67 


36 


0.205±0.105±0.006 


0.278 






s'/s > 0.7225 


42.5 


16.5 


0.440±0.138±0.011 


0.556 




"e~ 


1 cos 616- 1 < 0.7 
and ^acoi < 10° 


595 


69 


0.795±0.024±0.004 


0.815 






Asymmetries at -y/s 


- 205 GeV 






Nf 


iVb 




/iSM 

ID 




V- 


s'/s > 0.01 


270 


148 


0.234±0.051±0.006 


0.277 






s'/s > 0.7225 


160 


52 


0.512±0.063±0.005 


0.553 




s'/s > 0.01 


133 


58 


0.317±0.076±0.006 


0.277 






s'/s > 0.7225 


87 


26 


0.575±0.092±0.011 


0.553 




'e 


cos6'c-| < 0.7 
and 6'acoi < 10° 


1248 


147 


0.792±0.016±0.004 


0.816 






Asymmetries at y/s 


- 207 GeV 






N{ 


A^b 




/iSM 






s'/s > 0.01 


488 


251 


0.264±0.038±0.006 


0.277 






s'/s > 0.7225 


261 


86 


0.508±0.050±0.005 


0.552 


T" 




s'/s > 0.01 


205 


105 


0.230±0.061±0.007 


0.277 






s'/s > 0.7225 


142 


53 


0.472±0.075±0.011 


0.551 






1 cos6'c-| < 0.7 
and 6'acoi < 10° 


2062 


216 


0.814±0.012±0.004 


0.816 



Table 6: Continued 
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cos^l 



qq 

d(T/d| cos ^1 / pb 



189 GeV 



192 GeV 



196 GcV 



200 GcV 



202 GeV 



205 GcV 



207 GeV 



0.0, 
0.1, 
0.2, 
0.3, 
0.4, 
0.5, 
0.6, 
0.7, 
0.8, 
0.9, 



0.1] 
0.2] 
0.3] 
0.4] 
0.5] 
0.6] 
0.7] 
0.8] 
0.9] 
1.0] 



17.4±1.0 
17.9±1.0 
17.8±1.0 
18.3±1.1 
18.6±1.1 
21.3±1.1 
24.8±1.2 
25.8±1.2 
26.6±1.3 
31.0±1.7 



20.8±2.8 
14.6±2.4 
18.1±2.6 
17.6±2.6 
23.6±3.0 
23.2±3.0 
23.1±3.0 
23.8±3.0 
29.1±3.3 
26.2±4.0 



16.4±1.5 
14.4±1.5 
14.4±1.5 
16.1±1.6 
19.2±1.7 
20.5±1.7 
23.4±1.8 
21.3±1.7 
23.0±1.8 
28.9±2.6 



11.8±1.3 
15.8±1.5 
11.9±1.3 
16.2±1.5 
17.4±1.6 
18.0±1.6 
20.8±1.7 
21.4±1.7 
27.6±2.0 
27.7±2.5 



13.6±2.0 
16.1±2.2 
15.4±2.2 
20.5±2.5 
16.3±2.2 
14.4d=2.1 
18.6±2.4 
21.1±2.5 
23.4±2.6 
25.3±3.4 



12.3±1.3 
12.8±1.4 
12.6±1.4 
15.8±1.5 
16.5±1.5 
19.1±1.6 
20.4±1.7 
22.4±1.7 
21.7±1.7 
28.4±2.5 



13.0±1.0 
13.4±1.1 
10.7±1.0 
16.6±1.2 
15.0±1.1 
17.1±1.2 
18.9±1.3 
17.6±1.2 
20.9±1.3 
24.7±1.8 



Table 7: Differential cross-sections for qq production. The values are for s' / s > 0.7225 and are corrected to no 
initial- and final-state radiation. Errors are statistical only; systematic errors are given in Table 11. 



interference between 



COS 9 



dfj/dcos 9 / ph 



189 GeV 



192 GeV 



200 GeV 


202 GeV 


205 GeV 


207 GeV 


n oil 0.45 
U.Oi=Co.28 
Q 17-1-0.23 
="=0.13 

0.56±[JJ4 
0.60±0:i 
0.9±0.3 
1.3±0.3 


0.3±oj 


n 7«-|_0.58 
u. /uznQ 42 

Q 4Q-U0.30 

n 72-^0.35 

o.60±<]:2l 

^■'-'='=0.3 

1.5±0.3 


n 44-U0.34 

U. ^^3=0.25 

0.49±|]:22 
n 72-^0.26 

' ^="=0 20 
U. /UZCq 20 

1.01±0.20 

0.86±0.19 


1.9±0.4 




2.2±0.4 


1.7±0.3 


1.6±0.3 


-'-•^='=0.5 


1.8±0.4 


2.0±0.3 


2.8±0.5 


2-8±o:7 


2.2±0.4 


2.6±0.3 


3.8±0.6 




3.2±0.6 


3.3±0.5 



-1.0,-0.8] 
-0.8,-0.6] 
-0.6,-0.4] 
-0.4, -0.21 



-0.2, 
0.0, 
0.2, 
0.4, 
0.6, 
0.8, 



0.0] 

0.2] 
0.4] 
0.6] 
0.8] 
1.0] 



0.67±°:i 

0.47±°:1« 

0.54±0:?o 

0.62±0.13 

1.27±0.19 

1.35±0.20 

2.03±0.25 

2.15±0.26 

2.85±0.30 

3.77±0.42 



1 o I 1.3 
-'-•'J ="=0.8 

o.7±g:I 

0.5±oj 

0.8±n 
2.0±[]-? 
1.9±rr 

^- ' ="=0.8 
-'-•^='=0.9 



0.19+0;?^ 
0.28+J]:?^ 
0.79±{]:i 
0.54±S:i 
0.54±gi? 
1.7+0.3 
1.6+0.3 
2.2+0.4 
3.0+0.5 
3.9+0.6 



Table 8: Differential cross-sections for //"'"/x production. The values are for s'/s > 0.7225 and are corrected to no interference between 
initial- and final-state radiation. Errors are statistical only; systematic errors are given in Table 11. 



T+T 



COS^ 



do"/dcos 9 / ph 



189 GeV 



192 GeV 



189 GeV 200 GeV 202 GeV 



205 GeV 



207 GeV 



-1.0,-0.8] 
-0.8,-0.6] 
-0.6,-0.4] 
-0.4,-0.2] 



-0.2, 
0.0, 
0.2, 
0.4, 
0.6, 
0.8, 



0.0] 
0.2] 
0.4] 
0.6] 
0.8] 
1.0] 



-'-•-'-=1=0,8 

1.6d=0.3 
2.0±0.3 
2.7±0.4 
3.1±0.5 
4.7±0.8 



n 1 _|_1.6 
'-'•-'-='=0.0 

-0 l±°-5 
"-'•-'-='=0.0 

o.9±i:? 

3.7±}i 



0.7±n 
0.9±°i 

-'-•-'-='=0.5 

1 3+0-^ 

-'-•'J='=0.5 

2 0±°-S 

^•"-'='=0.6 

1.9±0.5 
4.1±?:0 



o.od=0:t 

0.5±oi 
0.8±oj 

2.1±0.5 
2.4±0.5 
3.4±0.7 



2.5±t 
0.3±l 
0.6±l 
1.2±l 

0.9±l 
1.2±l 



3.0±^ 



4.0 
2.6 


-0 1 +0-6 

'-'•-'-='=0.0 


1.0 
0.5 


7+°-6 

U. ( 11=0.4 


1.0 
0.6 




1.2 
0.8 


0.4±n 


1.0 
0.6 


'-'•^='=0.4 


1.1 
0.6 


-'-•^=1=0.5 


1.1 
0.7 


-'-•'-'='=0.6 


1.4 
1.1 


1 9±°-'^ 

-'-•^='=0.6 


1.4 
1.1 


2.5±0.6 


3.0 
2.0 


3.7±?:0 



o.2±n 

0.8±oj 

o.5±0:| 

-'-■^='=0.4 

1.2±0.3 
1.4±0.3 
2.6±0.4 
2.0±0.4 
2.5±\i 



Table 9: Differential cross-sections for r'^r production. The values are for s'/s > 0.7225 and are corrected to no 
initial- and final-state radiation. Errors are statistical only; systematic errors are given in Table 11. 



interference between 



e+e 


cos^ 


dfj/dcos^^ / pb 




189 GeV 


192 GeV 


196 GeV 


200 GeV 


202 GeV 


205 GeV 


207 GeV 


[-0.90, 


-0.72] 


1.3±0.2 


-'-•^='=0.5 


-'-•'J ="=0.3 


1.5±0.3 


^•^=■=0.6 


0.8±[ii 


1.4±0.2 


[-0.72, 


-0.54] 


2.1±0.3 


O C 1 Q 9 

2.5±q;7 


1.5±0.3 


1.8±0.4 


-'-•'J ="=0.4 


1.7ib0.4 


1.7±0.3 


[-0.54, 


-0.36] 


2.4±0.3 




2.2±0.4 


2.0±0.4 




2.3±0.4 


1.7±0.3 


[-0.36, 


-0.18] 


2.5±0.3 




2.9±0.5 


2.5±0.4 


3.8±0.8 


2.3±0.4 


2.8±0.3 


[-0.18, 


0.00] 


3.8±0.4 


^•^=■=0.8 


3.5±0.5 


3.9±0.5 




3.8±0.5 


2.8±0.4 


[ 0.00, 


0.09] 


4.5±0.5 




5.1±0.9 


5.1±0.9 


5.0±\i 


5.7±0.9 


3.9±0.6 


r n no 
[ u.uy. 


n 1 Ql 
U.ioJ 


6.3±0.6 


6.6±ll 


5.6±0.9 


6.2±1.0 


8.2±1.6 


6.9±1.0 


6.4±0.7 


[ 0.18, 


0.27] 


Q 9-l-n 8 


Q O-UI Q 


8 A -1-1 1 


Q 7-1-1 9 


o Q_|_i a 
o.yici.o 


7 94-1 n 


7 /i-Ln 8 


[ 0.27, 


0.36] 


12.3±0.9 


13.7±2.3 


12.7±1.4 


10.0±1.2 


10.8±1.8 


11.2±1.3 


11.9±1.0 


[ 0.36, 


0.45] 


19.7±1.1 


21.2±2.9 


14.9±1.5 


15.6±1.5 


14.6±2.1 


14.2±1.4 


17.1±1.2 


[ 0.45, 


0.54] 


31.9±1.4 


31.0±3.5 


26.7±2.0 


28.8±2.1 


24.5±2.7 


27.2±2.0 


25.1±1.5 


[ 0.54, 


0.63] 


51.6±1.8 


45.2±4.3 


50.1±2.8 


48.4±2.7 


45.4±3.8 


43.8±2.5 


41.9±1.9 


[ 0.63, 


0.72] 


94.9±2.4 


91.6±5.9 


90.5±3.6 


87.1±3.5 


84.4±5.0 


79.0±3.3 


77.3±2.5 


[ 0.72, 


0.81] 


215±4 


206±9 


199±6 


194±5 


184±8 


173±5 


174±4 


[ 0.81, 


0.90] 


684±7 


663±16 


640±10 


606±10 


606±14 


578±9 


565±7 



Table 10: Differential cross-sections for e"^e production for ^acoi < 10°. Errors are statistical only; systematic errors are given in 
Table 11. 



v/i / GcV 


189 


192 


196 


200 


202 


205 


207 


qq 


b 1 b ^ U.Ui 


0.90 


0.92 


0.97 


1.00 


1.00 


1.00 


1.04 






0.57 


0.54 


0.53 


0.52 


0.52 


0.51 


0.54 




b j b ^ U.UI 


1.12 


1.13 


1.13 


1.14 


1.19 


1.17 


1.16 




o / o ^ U. / ZZO 


1.08 


1.08 


1.07 


1.07 


1.08 


1.08 


1.08 




s'/s > 0.01 


2.6 


2.6 


2.6 


2.6 


2.5 


2.6 


2.6 




s' Is > 0.7225 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


e+e~ 


A: |cOS^e±| < 0.9, ^acol < 170° 


0.51 


0.46 


0.46 


0.46 


0.46 


0.59 


0.59 




B: 1 cos^c-l < 0.7, ^acoi < 10° 


0.31 


0.36 


0.36 


0.36 


0.36 


0.39 


0.39 




C: |cOS^e± < 0.96, ^acol < 10° 


0.44 


0.43 


0.43 


0.43 


0.43 


0.44 


0.44 




D: -0.9 < COS^e- < -0.7, ^acol < 10° 


1.10 


1.16 


1.16 


1.16 


1.16 


1.18 


1.18 




D: |cOS^e-| < 0.7, ^acol < 10° 


0.47 


0.50 


0.50 


0.50 


0.50 


0.52 


0.52 




D: 0.7 < cos^, < 0.9. 0,,,„i < 10° 


0.77 


0.86 


0.86 


0.86 


0.86 


0.89 


0.89 



Table 11: Total systematic errors, in %, excluding those on the luminosity measurement, for 
each channel at each nominal centre-of-mass energy. For the e''"e~ final state, the first three 
rows refer to the total cross-section measurements, while the last three refer to the differential 
cross-section, for which extra cuts are applied. For the hadronic, and t~^t~ final states, 

the values given for s'/s > 0.7225 apply to both the total cross-section and the differential 
cross-section. 



qq 




s'/s > 0.01 


s'/s > 0.7225 


MC statistics (efficiency) 


0.03 


0.06 


MC statistics (background) 


0.03 


0.04 


ISR modelling 


0.12 


0.02 


Fragmentation modelling 


0.37 


0.26 


Detector effects 


0.19 


0.15 


s' determination 


0.03 


0.19 


W"'"W^ rejection cuts 


0.07 


0.13 


W+W" background 


0.21 


0.26 


Other background 


0.87 


0.19 


Interference 


0.04 


0.13 


Total 


1.00 


0.52 



Table 12: Systematic errors, in %, on the hadronic cross-section measurements at 200 GeV. 
Values at other energies are very similar, the total errors are given in Table 11. Errors on the 
luminosity measurement are given in Table 3. 
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s'/s > 0.01 


s'/s > 0.7225 


ML statistics (emciencyj 


0.08 


0.08 


MC statistics (background) 


0.10 


0.08 


MC statistics (feedthrough) 




0.03 


Efficiency 


1.00 


1.00 


ISR modelling 


0.20 


0.07 


Feedthrough 




0.16 


Cosmic background 


0.20 


0.20 


Other background 


0.46 


0.20 


Interference 


0.01 


0.12 


Total 


1.11 


1.07 



Table 13: Systematic errors, in %, on the cross-section measurements at 200 GeV. The 

errors at other energies are very similar, the totals at each energy are given in Table 11. Errors 
on the luminosity measurement are given in Table 3. 



r 








s'/s > 0.01 


s'/s > 0.7225 


MC statistics (efficiency) 


0.19 


0.21 


MC statistics (background) 


0.22 


0.20 


MC statistics (feedthrough) 




0.05 


Efficiency 


2.34 


2.34 


ISR 


0.20 


0.04 


Feedthrough 




0.15 


Background 


0.99 


0.69 


Interference 


0.04 


0.01 


Total 


2.6 


2.5 



Table 14: Systematic errors, in %, on the r+r" cross-section measurements at 200 GeV. The 
errors at other energies are very similar, the totals at each energy are given in Table 11. Errors 
on the luminosity measurement are given in Table 3. 
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e+e 




A 


B 


c 




1 cos 6*6*1 < 0.9 


1 cos 6*6-1 < 0.7 


1 COS 6*6*1 < 0.96 




^acol < 170° 


^acol < 10° 


^acol < 10° 


MC statistics 


0.02 


0.04 


0.02 


4-fermion correction 


0.06 


- 


- 


Multiplicity cuts 


0.09 


0.04 


0.03 


Calorimeter energy scale/resolution 


0.01 


< 0.01 


0.08 


Two track requirement 


0.32 


0.30 


- 


Acceptance 


0.19 


0.14 


0.39 


Background 


0.25 


0.13 


0.17 


Total 


0. 16 


0.36 


0.13 




D 


D 


D 




-0.9 < cos6ie- < -0.7 


|cos6ie-| < 0.7 


+0.7 < cos0^^ < +0.9 




^acol < 10° 


6'acol < 10° 


^acol < 10° 


Multiplicity cuts 


0.05 


0.04 


0.05 


Calorimeter energy scale/resolution 


0.01 


0.01 


0.01 


Two track requirement 


0.48 


0.30 


0.48 


Opposite charge requirement 


0.64 


0.37 


0.64 


Charge misassignment 


0.50 






Acceptance 


0.30 


0.10 


0.30 


Background 


0.60 


0.11 


0.02 


Total 


1.16 


0.50 


0.86 



Table 15: Systematic errors, in %, on the e"'"e~ cross-section and angular distribution measure- 
ments at 200 GeV. Values at other energies are very similar, the totals at each energy are given 
in Table 11. Errors on the luminosity measurement are given in Table 3. In the case of the 
angular distribution, acceptance D, the errors arising from Monte Carlo statistics are included 
in the statistical errors given in Table 10. 
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otanaara Moael 


Data / Standard Model 






-.,2 / J ^ f 

X / d.o.t 


Mean 


^ 1 d.o.f 




s / s > U.Ul 


5.9/7 


0.997±0.010 


5.8/6 




/ ^ ^ n '700 

s /s > OJZZb 


5. (J/7 


0.990±0.011 


4.2/6 


^( ,.+ ..— ^ 


s /s > U.Ui 


9. 8/7 


1.012±0.021 


9.5/6 




/ / 700 r 

s /s > U.7225 


/I 7 /7 

4.7/7 


U.994±U.U28 


4.6/6 


A 1 + — \ 


s /s > U.Ul 


II.d/7 


U.975±U.U65 


11.5/6 




s /s > [j.izZo 


lz.d/7 


U.999±U.U4U 


12.3/6 




s fs > U.Ul 


z.d/7 


l.U45±U.U39 


1.3/6 




s /s > 0.7225 


2A/1 


l.U52±U.U44 


1.1/6 


/I /"_+_— ^ 
Aq,[t^t ) 


s /s > U.Ul 


P O /7 

o.o/7 


i.Uio±U.U9o 


o.z/o 




s'/s > n 7225 


13 8/7 


l.U33±U.U57 


13.5/6 


a {£+£-) 


s'/s > U.Ul 


12.4/14 


l.U2U±U.U19 


11.3/13 




s'/s > u.7225 


7.1/14 


1.U1U±U.U24 


6.9/13 




s'/s > U.Ul 


17.9/14 


U.988±U.U54 


17.9/13 




,s7,s > 0.7225 


2(3.2/11 


1.010±0.033 


26.1/13 



Table 16: Comparison of measurements with Standard Model predictions. The first column 
gives the value of the measured cross-sections or asymmetry values at 189-207 GeV presented 
here with respect to the Standard Model predictions. The second and third columns give the 
results of fits to the mean ratios of data to Standard Model predictions. Values for t^l" are 
for and r+r" together. 
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Using cross-sections 






Using ratios 






Fit 


Standard Model 


Fit 


Standard Model 


/ GeV 


l/'-^cm 


xVd-o.f. 






l/'-*^cm 






xVd-o.f. 


188.6 


127.1±|J 


2.0/4 


127.9 


2.0/5 


126.2±|:2 


1.7/3 


127.9 


1.9/4 


191.6 


134.9±l^:f 


10.0/4 


127.9 


10.4/5 


138.1±\t,i 


9.1/3 


127.9 


10.0/4 


195.5 


131.2±i;? 


8.4/4 


127.9 


8.7/5 


129.9±^i 


8.1/3 


127.9 


8.2/4 


199.5 


130.3±^i 


6.4/4 


127.8 


6.6/5 


130.3±|° 


6.4/3 


127.8 


6.5/4 


201.6 


134.1±^^5' 


2.8/4 


127.8 


3.3/5 




2.8/3 


127.8 


3.1/4 


204.9 


122.1±|5 


0.7/4 


127.8 


1.8/5 


123.7±|| 


0.6/3 


127.8 


0.9/4 


206.6 


123.4±|o 


4.6/4 


127.8 


5.7/5 


117.9±t^ 


0.6/3 


127.8 


4.5/4 


193.2 


127.4±i:J 


59.1/59 


127.9 


59.1/60 


126.7±i:| 


50.4/47 


127.9 


50.7/48 



Table 17: Results of fits for ctem- The first seven rows show the fits to data at each energy, the last row the combined fit to these 
data and measurements at 130-183 GeV [2,3]. The Standard Model values of l/cicm, and the between the measurements and the 
Standard Model predictions are also given for comparison. Results are shown for the fits using cross-sections, and also for the fits to 
cross-section ratios, as discussed in the text. 



Parameter 


Without lepton 
universality 


With lepton 
universality 


Standard Model 
prediction 


mz / GeV 


91.1882±0.0033 


91.1872±0.0033 




Tz / GeV 


2.4945±0.0041 


2.4943±0.0041 


r> Ar\r<r\ +0 0(116 

^■^^"^-0.0029 


^tot 
'had 

J had 


2.963±0.009 
0.131±0.078 


2.963±0.009 
0.144±0.078 


_ _ _ ^ _i_n 00*^7 

2-9650+;:S 


' e 


0.14134±0. 00069 






y. LU L 


0.14215±0. 00056 






^tot 

T 


0.14228±0. 00074 






„tot 




0.14199±0.00050 




„tot 
Je 


-0.080±0.044 






tot 
J fi 


-0.008±0.019 






^'tot 
Jt 


-0.004±0.025 






•tot 




-0.014±0.015 


0.00439+JX° 


e 


0.00138±0. 00084 








0.00270±0. 00043 






„fb 
' r 


0.00248±0.00057 










0.00243±0.00032 


0.00280 


■fb 
Je 


U. /00±U.U lO 






.•fb 
J IX 


0.758±0.024 






jib 
Jt 


0.788±0.030 










0.767±0.018 


7087+0-0005 

u. < »0( _Q 0006 


xVd.o.f. 


207.0 / 247 


213.3 / 255 





Table 18: Results from the fit to all LEP 1 and LEP2 data for the S-matrix parameters with 
and without the assumption of lepton universality. The last column gives the Standard Model 
predictions. An S-matrix fit to only LEP 1 data without lepton universality gives a of 146.6 
with 187 d.o.f. [39]. 
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Parameter 


1 2 


3 4 


5 6 7 


8 9 10 


11 12 13 


14 15 16 


1 mz 

2 Tz 


1.00 0.04 
0.04 1.00 


0.06 -0.40 
0.92 -0.08 


-0.08 0.04 0.03 
0.57 0.71 0.54 


-0.11 -0.12 -0.11 
-0.03 0.01 0.02 


-0.06 0.07 0.06 
0.00 0.02 0.01 


0.01 -0.04 -0.03 
0.00 0.04 0.04 


q „tot 
"J '^had 

4 J^d 


0.06 0.92 
-0.40 -0.08 


1.00 -0.09 
-0.09 1.00 


0.57 0.71 0.54 
0.00 -0.07 -0.05 


-0.04 0.00 0.01 
0.08 0.09 0.08 


0.01 0.02 0.02 
0.03 -0.04 -0.03 


0.00 0.04 0.04 
0.00 0.03 0.02 


5 r*°* 

c „tot 
^^ 


-0.08 0.57 
0.04 0.71 
0.03 0.54 


0.57 0.00 
0.71 -0.07 
0.54 -0.05 


1.00 0.45 0.33 
0.45 1.00 0.41 
0.33 0.41 1.00 


0.03 0.02 0.02 
-0.03 0.08 0.01 
-0.02 0.00 0.09 


0.17 -0.01 -0.01 
0.01 0.03 0.01 
0.00 0.01 0.03 


-0.02 0.03 0.03 
0.00 0.09 0.03 
0.00 0.03 0.10 


8 jf* 

9 jT 
10 


-0.11 -0.03 
-0.12 0.01 
-0.11 0.02 


-0.04 0.08 
0.00 0.09 
0.01 0.08 


0.03 -0.03 -0.02 
0.02 0.08 0.00 
0.02 0.01 0.09 


1.00 0.02 0.02 
0.02 1.00 0.03 
0.02 0.03 1.00 


0.00 -0.01 -0.01 
0.01 0.06 -0.01 
0.01 -0.01 0.06 


0.24 0.01 0.00 
0.00 0.33 0.01 
0.00 0.01 0.29 


11 
12 

13 rf 


-0.06 0.00 
0.07 0.02 
0.06 0.01 


0.01 0.03 
0.02 -0.04 
0.02 -0.03 


0.17 0.01 0.00 
-0.01 0.03 0.01 
-0.01 0.01 0.03 


0.00 0.01 0.01 
-0.01 0.06 -0.01 
-0.01 -0.01 0.06 


1.00 -0.01 -0.01 
-0.01 1.00 0.02 
-0.01 0.02 1.00 


0.06 0.00 0.00 
0.00 0.12 0.00 
0.00 0.00 0.12 


14 

15 J? 
16 


0.01 0.00 
-0.04 0.04 
-0.03 0.04 


0.00 0.00 
0.04 0.03 
0.04 0.02 


-0.02 0.00 0.00 
0.03 0.09 0.03 
0.03 0.03 0.10 


0.24 0.00 0.00 
0.01 0.33 0.01 
0.00 0.01 0.29 


0.06 0.00 0.00 
0.00 0.12 0.00 
0.00 0.00 0.12 


1.00 0.00 0.00 
0.00 1.00 0.02 
0.00 0.02 1.00 



Table 19: Error correlation matrix for the S-matrix fit without lepton universality. 



Parameter 



6 



8 



1 
2 

3 
4 
5 



^tot 
'had 

■tot 
J had 

„tot 



6 fe 



tot 



.fb 



8 



fb 



1.00 
0.04 
0.06 

-0.39 
0.01 

-0.19 
0.07 

-0.06 



0.04 
1.00 
0.92 
-0.08 
0.80 
0.00 
0.02 
0.06 



0.06 
0.92 
1.00 

-0.09 
0.80 

-0.01 
0.03 
0.05 



-0.39 
-0.08 
-0.09 

1.00 
-0.06 

0.14 
-0.03 

0.04 



0.01 
0.80 
0.80 
-0.06 
1.00 
0.06 
0.05 
0.09 



-0.19 
0.00 

-0.01 
0.14 
0.06 
1.00 
0.05 
0.31 



0.07 
0.02 
0.03 
-0.03 
0.05 
0.05 
1.00 
0.11 



-0.06 
0.06 
0.05 
0.04 
0.09 
0.31 
0.11 
1.00 



Table 20: Error correlation matrix for the S- matrix fit assuming lepton universality. 
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Channel 




LL 


RR 


LR 


RL 


W 


AA 


LL+RR 


LR+RL 


Cdb 






[±1,0,0,0] 


[0, ±1,0,0] 


[0,0, ±1,0] 


[0,0, 0,±1] 


[±i,±i,±i,±i] 


[±l,±l,Tl,Tl] 


[±1, ±1,0,0] 


[0,0,±1,±1] 


[±i,±i,±l,±i] 




A+ 


— 0.017 

4.7 


009+° °}° 
4.7 


-0 009+" °J,2 

u.uu.y_0.009 
8.1 


-0 009+" °J,S 

" ""'^-O.OOO 

8.1 


-0 002+"-°"} 

" -0.004 

12.6 


006+° °"-^ 

" """-0.005 

8.1 


004+!!-°"° 

"""— 0.008 

6.8 


-0 005+°-!;°-^ 

" """-0.005 

11.8 


-0 003+!!-°"^ 

" """-0.007 

9.2 




A_ 


6.1 


6.0 


6.2 


6.2 


10.6 


11.9 


8.5 


8.6 


7.9 




En 


-0.002+°°°^ 
8.1 


-0 002+" °"° 

"•""^ — 0.009 

7.7 


002+° "}2 
6.3 


002+° °}2 

"•""^ — 0.013 

6.3 


000+° °°-^ 

— 0.003 

12.7 


-0 001+° °°! 

u.uvj-_o.004 

11.3 


-0 ooi+°-!!°i; 

"•""-^—0.004 

11.1 


001+° °"° 

"•""^ — 0.006 

8.7 


000+° °°° 

"•""" — 0.006 

9.6 




A_ 


7.3 


7.0 


6.3 


6.3 


12.4 


10.2 


10.1 


9.3 


9.2 




Fn 

A+ 


"•'-'-^"—0.014 

4.9 


017+° °" 

u.ui 1 —0.015 

4.7 


-0 002+" °3,^ 

"•""^—0.020 

5.7 


-0 002+°-°i^ 

"■""^—0.020 

5.7 


004+° °°^ 
8.6 


008+!l!!i!l 

0.007 

6.7 


008+!!°°l 

0.007 

6.8 


-0 ooi+!i;!V^ 

"•""-^—0.010 

7.8 


007+!!!!!!?: 

"•"" ' —0.009 

6.4 




A_ 


7.2 


6.9 


4.6 


4.6 


11.1 


10.7 


10.0 


7.0 


8.4 


i+i- 


A+ 


0.004+°-°°? 

'-'•'-"-'^—0.007 

7.7 


0.004+°°°? 

'-'•'-"-'^—0.007 

7.4 


-0.005+S-S^I 
9.3 


-0.005+^-S^I 
9.3 


0.000+S-^S? 
15.2 


0.003+S-^S? 

v^.ww(j_0 003 

10.5 


0.002+°-°°l! 

w.v^i^^_0 004 

10.6 


-0.002+°°°! 

v^.ww^_0 004 

13.2 


0.000+S-^R! 

i^.www_o 004 

11.2 




A_ 


9.5 


9.2 


7.3 


7.3 


15.1 


15.4 


13.3 


10.3 


11.4 




A+ 


-0 021+° °2i 
8.2 


0.018+°°23 
4.3 


0.007+°"?" 
4.9 


0.011+°°}^ 

— 0.010 

3.1 


0.011+°°!? 

— 0.012 

5.7 


-0.009+°°°° 

\j.vjvj,y_0 025 

12.0 


-0.007+°°!! 

7.5 


010+° °^^ 

vj.»j-nj_o 009 

3.9 


0.011+°°?? 

— 0.010 

3.7 




A_ 


3.7 


7.0 


6.1 


9.3 


10.4 


5.0 


5.7 


10.3 


9.7 


combined. 


A+ 


ooi+l!-^;}^ 

U.WU J._Q QQg 

9.2 


'^•'-"-"-'—0.007 

7.2 


-0 oo3+ii;;iii 

9.4 


ooi+^ji;}^ 

9.0 


OOQ+anai 
u.wuu_o,o02 

14.7 


001+!! ;!1!2 

12.6 


002+° °"^ 

w.uu^_0_003 

11.1 


-0 ooi+!!!!;!2 

u.ww-L_o_o03 

12.9 


002+° °°! 

u.ww^_0_004 

10.6 




A_ 


9.4 


10.1 


7.8 


10.1 


16.2 


14.9 


13.2 


12.3 


13.0 


uu 


A+ 


0.009+°°°^ 

"•"""^ — 0.009 

5.9 


014+"-"}-^ 

"•"^^ — 0.013 

4.8 


048+°-°!^ 

"•"^" — 0.047 

2.4 


018+° °S 

"•"-^" — 0.044 

3.3 


005+" "°^ 

— 0.004 

8.5 


007+°-!!°? 

' —0.007 

6.7 


0.006+°°"t 
7.8 


0.038+°°^? 

w.»jtj"j_o 038 

3.0 


009+° °°° 

"•"""^ — 0.009 

6.0 




A_ 


9.1 


7.7 


5.4 


4.1 


13.0 


10.7 


11.8 


6.2 


9.3 


dd 


A+ 


-0 011+"°!! 

"•"^^-0.011 

8.6 


-0 033+" ?-^^ 
6.0 


-0 036+" °^;i 

v.uou_o.049 

4.6 


045+° °?? 

"•"^^ — 0.043 

2.4 


-0 009+° °?° 

"•""-^—0.010 

9.9 


-0 007+° °°! 

' —0.008 

10.9 


-0 008+"-!!°^ 

"•"""-0.008 

10.2 


019+° °?? 

"•"^■^-0.032 

3.7 


-0 034+° "3? 

"•""^-0.085 

6.3 




A_ 


5.5 


2.2 


2.9 


5.4 


5.6 


6.7 


6.5 


5.2 


2.7 


uu + dd 


So 
A+ 


0.0181°;°- 
3.7 


3.0 


0.0181°;°- 
3.7 


3.0 


01 1 +° "^° 

U.Uii_o.oil 

3.8 


0.0001°;°- 
6.1 


0.028l°;°or7 
3.6 


U.UZ/_o 027 

3.6 


n ni R+0.025 
O.OI6I0.015 

2.8 




A_ 


5.1 


6.2 


5.2 


6.2 


9.7 


6.1 


7.3 


7.3 


7.9 



Table 21: Results of the contact interaction fits to the non-radiative hadron and lepton-pair data. The numbers in square brackets are 
the values of [?7ll,''7rr,''7lr, ''7rl] which define the models. £0 is the fitted value of £ = 1/A^, A-t are the 95% confidence level limits; the 
values for A+ and A_ correspond to the upper and lower signs, respectively, of the rjij values. The units of A are TeV, those of £0 are 
TeV~l 



Model: 


X V' ?7 LR SSM 


m'zT / GeV 
O'^ / mrad 
1 mrad 


781 366 515 518 1018 
1.94 2.58 3.31 1.90 0.91 
-0.99 -1.29 -4.47 -0.98 -4.22 



Table 22: Onc-dimensional limits at 95% confidence level on the Z' mass, m^™, and the mixing 
angle, and ^ for various Z' models. The Z mass is free during the fit and the other three 
Standard Model parameters (ttg, ''^top and mniggs) are fixed at their default values. 
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Figure 1: Integrated luminosity collected by OPAL, and used in these analyses, during 1998, 
1999 and 2000. The dashed lines indicate the division of the 2000 data into the two centre-of- 
mass energy bins, 202.5 GeV< ^/s <205.5 GeV and ^/s >205.5 GeV. The precise amount of 
data used in each analysis varies slightly from channel to channel. 
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Figure 2: (a) The distribution of the difference in radial coordinate between the two clusters in 
Bhabha scattering events used for the silicon-tungsten luminosity measurement. Distributions 
of the radial coordinates of clusters are shown for (b) the 'narrow' side and (c) the 'wide' side 
calorimeter. Distributions arc shown after all cuts except the acollinearity cut in (a) and the 
inner and outer radial acceptance cuts, on that side, in (b). Points show the data taken in the 
year 2000, while the histograms show the Monte Carlo expectation. The vertical bars show the 
positions of the cuts which define the acceptance, with the arrows pointing into the accepted 
region. 
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Figure 3: The pad boundary image at the inner acceptance cut (i?in = 9.45 cm) used in the 
luminosity measurement. The nominal pad boundary is conventionally set at zero. The points 
show the fraction of events with pad maximum beyond the nominal cut as a function of distance 
from the pad boundary for the layer located after 7 radiation lengths, for data taken in 2000. 
The solid curves show the fitted functions used to determine the coordinate offsets. 



60 




Figure 4: Distributions of variables used in the selection of hadronic events: (a) number of 
tracks, (b) number of electromagnetic calorimeter clusters, (c) ratio of the visible energy to the 
centre-of-mass energy and (d) energy balance along the beam direction. The points show the 
data for all centre-of mass energies combined and the histograms the Monte Carlo predictions 
normalized to the integrated luminosity of the data. In each case the distribution is shown after 
the selection cuts associated with the other three variables have been applied. The positions of 
these cuts are indicated by the vertical bars, with the arrow pointing into the accepted region. 
The W"^W~ rejection cuts have not been applied. Background labelled 'other' is mainly r+r^. 
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Figure 5: Distributions of the QCD matrix element for four-jet production W420 for events 
(a) passing the hadronic event selection and the non-radiative s' cut before applying the 
W+W" veto, (b) additionally failing the W"'"W~ qqqq veto or (c) additionally passing 
the W^W~ qqqq veto. The points show the data for all centre-of-mass energies combined 
and the histograms the Monte Carlo predictions normalized to the integrated luminosity of the 
data. The contributions from signal events and from the various sources of background are 
indicated, while the fit region (discussed in the text) is shown by the arrows. 
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Figure 6: (a) Time difference between baclc-to-baclc hits in the time-of-ffight counters in the 
barrel region. Events in the combined data sample which pass all selection criteria 

except for the cosmic veto are included, if they have back-to-back TOF hits, (b) Ratio of the 
visible energy, defined as the sum of the muon momenta plus the energy of the highest energy 
electromagnetic calorimeter cluster, to the centre-of-mass energy, for //"'"/i" candidates passing 
all cuts except those on the visible energy and the mass of the muon pair. The points show 
the combined data and the histograms show the Monte Carlo expectation, normalized to the 
integrated luminosity of the data, with the background contributions as indicated. The vertical 
bars indicate the positions of the cuts (for a centre-of-mass energy of 200 GeV in (b)), with the 
arrow pointing into the accepted region in each case. 
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Figure 7: (a) Ratio of the total event energy, defined as the scalar sum of track momenta plus the 
energy of electromagnetic calorimeter clusters, to the centre-of-mass energy, for t~^t~ candidates 
passing all cuts except those on the visible energy and its track and cluster components, (b) 
The missing momentum, calculated using electromagnetic calorimeter clusters, divided by the 
centre-of-mass energy for t'^t~ events passing all cuts except those on the missing momentum 
and the cosine of the polar angle of its direction. In each case, the points show the combined 
data and the histograms show the Monte Carlo expectations, normalized to the integrated 
luminosity of the data, with the background contributions as indicated. The vertical bars 
indicate the positions of the cuts, with the arrow pointing into the accepted region in each 
case. Note that, in the case of the total event energy in (a), further cuts are placed on the 
separate track and cluster components at both low and high values. 
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Figure 8: (a) The acollinearity angle distribution for events satisfying the inclusive e"'"e~ — > e^e~ 
selection, in the acceptance region A, | cos6'e±| < 0.9. (b) The distribution of the ratio of total 
electromagnetic calorimeter energy to the centre-of-mass energy for e"'"e~ e"'"e~ events in 
acceptance region B, | cos6'e-| < 0.7 and 6'acoi < 10°. (c) The same distribution for the large 
acceptance region, C, | cos 6^6*1 < 0.96 and ^acoi < 10°. Distributions are shown after all cuts 
except the one on the variable plotted. In each case, the points show the combined data and 
the histograms the Monte Carlo expectations, normalized to the integrated luminosity of the 
data, with the background contributions as indicated. The vertical bars indicate the positions 
of the cuts in the displayed variable, with the arrow pointing into the accepted region in each 
case. 
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Figure 9: The distributions of reconstructed y s' j s for (a) hadronic events, (b) events, (c) 

T~^T~ events and (d) e"'"e~ events with | cos^^c±| < 0.9 and 6'acoi < 170°, for all data combined. In 
each case, the points show the data and the histogram the Monte Carlo prediction, normalized to 
the integrated luminosity of the data, with the contribution from events with true s'/ s > 0.7225 
shaded in (a), (b) and (c), and the contribution from events with 6'acoi < 10° shaded in (d). 
The vertical bars in (a), (b) and (c) show the position of the cut used to select 'non-radiative' 
events. 



66 



OPAL e"^e"^ hadrons 




200 

Vs / GeV 



Figure 10: Measured total cross-sections {s'/s > 0.01) for hadronic events at lower energies [2, 
3,39] and this analysis. Cross-section measurements for s'/s > 0.7225 from this analysis 
and from [2, 3] are also shown; the values at 161 GeV and 172 GeV have been corrected 
from s'/s > 0.8 to s'/s > 0.7225 by adding the prediction of Zfitter for this difference 
before plotting. The curves show the predictions of Zfitter. The insets show the percentage 
differences between the measured values and the Zfitter predictions for the high energy points. 
The error bars on the differences represent statistical errors only; the size of the experimental 
systematic error is indicated by the shaded band. 
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Figure 11: Measured total cross-sections (s'/s > 0.01) for events at lower energies [2,3,39] 

and this analysis. Cross-section measurements for s' j s > 0.7225 from this analysis and from [2, 
3] are also shown; the values at 161 GeV and 172 GeV have been corrected from s'/s > 0.8 
to s'/s > 0.7225 by adding the prediction of Zfitter for this difference before plotting. The 
curves show the predictions of Zfitter. The insets show the percentage differences between 
the measured values and the Zfitter predictions for the high energy points. The error bars 
on the differences represent statistical errors only; the size of the experimental systematic error 
is indicated by the shaded band. 
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Figure 12: Measured total cross-sections (s'/s > 0.01) for r+r^ events at lower energies [2,3,39], 
and this analysis. Cross-section measurements for s'/s > 0.7225 from this analysis and from [2, 
3] are also shown; the values at 161 GeV and 172 GeV have been corrected from s'/s > 0.8 
to s' I s > 0.7225 by adding the prediction of Zfitter for this difference before plotting. The 
curves show the predictions of Zfitter. The insets show the percentage differences between 
the measured values and the Zfitter predictions for the high energy points. The error bars 
on the differences represent statistical errors only; the size of the experimental systematic error 
is indicated by the shaded band. 
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Figure 13: Measured cross-sections for e"'"e~ events at lower energies [2,3,39], and this analysis. 
The curves show the predictions of Bhwide. The insets show the percentage differences between 
the measured values and the Bhwide predictions for the high energy points. The error bars on 
the differences represent statistical errors only; the size of the experimental systematic error is 
indicated by the shaded band. 
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Figure 14: (a) and (b) Measured asymmetries for inclusive {s'/s > 0.01) and non-radiative 
{s'/s > 0.7225) samples as functions of ^/s for /j^fi~ and t~^t~ events. The curves show 
Zfitter predictions for s'/s > 0.01 (solid) and s'/s > 0.7225 (dashed), (c) Measured forward- 
backward asymmetry for e+e~ with | cos 6*6-1 < 0.7 and 6'acoi < 10°, as a function of y/s. The 
curve shows the prediction of Bhwide. Lower energy data values are taken from [2,3, 39] for 
all channels. 
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Figure 15: Measured differential cross-sections for (a) hadronic events with s' j s > 0.7225 and (b) fi'^fi^ events with s'/s > 0.7225. The 
points show the luminosity-weighted average of all data from 189 GeV to 207 GeV, corrected to no interference between initial- and 
final-state radiation. The curves show the predictions of Zfitter without interference between initial- and final-state radiation (solid) 
and with interference (dashed). In each case the lower plot shows the ratio of the measurements to the Standard Model predictions 
(excluding interference) . 
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Figure 16: Measured differential cross-sections for (a) t~^t^ events witfi s'/s > 0.7225 and (b) e+e^ events witfi ^^acoi < 10°. Tlie points 
sliow tlie luminosity-weiglited average of all data from 189 GeV to 207 GeV, corrected to no interference between initial- and final-state 
radiation in (a). The curves in (a) show the predictions of Zfitter without interference between initial- and final-state radiation 
(sohd) and with interference (dashed). The curve in (b) shows the prediction of Bhwide. In each case the lower plot shows the ratio 
of the measurements to the Standard Model predictions (excluding interference in (a)). 
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Figure 17: Fitted values of l/aem as a function of Q, which is ^/s for the OPAL fits. The left plot shows the results of fits to OPAL 
data at each centre-of-mass energy and of the combined fit in which a^m runs with a slope obtained from fixing 1/Q;em(0) = 137.036. 
The right plot compares the results of the OPAL combined fits with values obtained by the TOPAZ experiment [43] and from fits 
to measurements of leptonic cross-sections and asymmetries at the DORIS, PEP, PETRA and TRISTAN c+c^ storage rings [44]. 
Measurements shown by open symbols rely on assuming the Standard Model running of ctcm for Qiumi below 4 GeV, whereas closed 
symbols indicate values derived from cross-section ratios which do not depend on luminosity, as discussed in Section 5.1. The solid line 
shows the Standard Model expectation, with the thickness representing the uncertainty, while the value of 1/Q;em(0) is shown by the 
dashed fine. 
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Figure 18: Confidence level contours in the mz - plane from the S-matrix fits with lepton 
universality. The dashed curve shows the 68% confidence level contour from the fit to LEP 1 
data alone, while the full and dotted curves show the 68% and 39% confidence level contours, 
respectively, from the fit to LEP 1 and LEP 2 data. The horizontal band indicates the Standard 
Model value of j^^. The vertical band is the la error on the Z mass from the five parameter 
fit [39] which should be compared with the 39% confidence level contour from the S-matrix fit. 
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Figure 19: 95% confidence level limits on the energy scale A resulting from the contact in- 
teraction fits to hadron and lepton-pair data. For each channel, the bars from top to bottom 
indicate the results for models LL to Obb in the order given in the key. The values for A+ and 
A_ correspond to the upper and lower signs, respectively, of the rjij values which define the 
models as given in Table 21. 
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Figure 20: Exclusion contours in the Z' mass - mixing angle plane at 95% confidence level 
for four Z' models. The Z mass is free during the fit and the other three Standard Model 
parameters {as, rritop and mniggs) are fixed at their default values. Leaving mtop and ctg free in 
the fit would lead to an increase of the width in by less than 10 %. 
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Figure 21: One dimensional limits at 95% confidence level for E6 and LR models as a function 
of tiie model angle. The particular cases of the rj, x, and LR symmetric models are indicated 
by the dots, (a) shows the limits on the Z' mass, and (b) shows the upper and lower limits 
on the mixing angle. They are obtained from a fit with a^, rritop and mniggs fixed but mz is 
free, (c) shows the absolute change in the limit on the mixing angle if as and mtop are free 
parameters but constrained by their experimental uncertainties, (d) shows the absolute change 
in the limit on the mixing angle if mniggs = 250 GeV instead of the default value of 115 GeV. 
In (c) and (d) the dashed curve denotes the change in the positive limit and the dotted curve 
denotes the change in the negative limit. 
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Figure 22: 95% confidence level exclusion contours on the axial and vector couplings of a Z' to 
leptons, for three values of the Z' mass. 
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